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Chapter I 


liTTH^C-nOS 

1-1 L'inear BonferronC Bound Probim 

The problem of estimating the probabilities of unions 
and intersections of events when all the information required 
about the event is not knovai, is one of the oldest unsolved 
pJ^oblems in probability theor^r. Generally the probability of | 

the objective event is to be estimated in the sense of giving ; 

its sharpest upper and lower bounds in terms of the avadlable i 

data. The linear Bonf err oni Bound problem is a problem of this [ 

hind. The origin of such type of problem goes as far back as | 

to Boole j~ • Sut it is only after the Vi?ork by Prechet CiiH. 1 
the problem has acquired its present status . Prechet de|ermines 
the sharpest upper and lovjer bounds for the proba-bilities of > 

unions and intersections in terms of their respective probabilities; 

1 

In view of the fact that Pre chefs result cannot be i 

used for estimating the probabilities of all compound events, 

several workers <[153, [l83> C 10 U > ClSl'l - C 19 H > 

11203. C293 - C313. Cssa. C343, C373> [1103, ' 

1347 * 3 , [3 5S3)> have been mot ivated, to extend the 

result to voider classes of problems. Linear Bonferroni Bound 

, i 

problem is also the result of such extensions. 
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P . • 
X 


Let A^, Ag,...,A^ Le any sequence of events. Let 
. = P(A. Oa. n ••• ) for all 1 < < ig < 

l^g***^v ^2 V 


... < < n. Consider the follov/ing linear inequality, 


5c + Z x-P^ + Z 

1 < i < n 1 < i^ < ig < n -^1-2 ^1^2 




1 < i^ < ig < ... i^ < n H^2 *’*^v ^ 1 ^ 2 ** v 

n 


<P(U A.) 
i=l ^ 


< x„ + Z x^- P - + Z X- ^ P.- - + ... 

" 1 < i < n ^ 1 £ ij_ < ig < n “^12 12 


Z X- • . p. . . , 

1 ^ i^ < ig < , . . < iy < n "^1^2 * ' ‘ ^1^2 * * ' \ 


( 1 . 1 . 1 ) 

Ihen the linear Bonferroni hound problem of degree v 


is to determine the real numbers x^, x^, x^^, x^, 


• • J -iX 


l-il 


r ^ ^2 * * * 1^ 


1^2' 

, such that the inequality (1.1,1) provides the 

n 


• r. 


i=l 


n 


the sharpest possible bounds of P( lli A-) in the similar manner 

i=l ■ 

is also called the linear Bonferroni bound problem. 

Though a classical problem in probability theory, there 
is no satisfactory solution of the linear Bonferroni bound 
problem up to now,^ At present the. exact computation of the 
bounds of degree two seems to be very difficult problem except 


sharpest possible bounds of P( <4 problem of determining 

. r ■ -i i , ■ ' ■ ■ ' ' 



3 


for fe?/ special cases, such as when P - = p -n t, 

1 1 ’ = ^ 12 » 

^ 1 i ^ 1 “-j ■fche bounds are readily avail-hip v, • 

j c.vaij_^,oie, liius searching 

for a sufficiently precise bound is still a phnii^v, • 

^ challenging problem. 


In the present work we restrict our s 


cope to finding the 


Bonferroni lower bound of degree two for Pf ^ 

■^vu A^). In the 

following section we give a brief account of'thp o • n -ui n 

01 tne available lower 

bounds for this problem. 

1.2 Knoim Bounds 

The most commonly knovm lower bound is 
n n 

P(_U A.) > 2 p _ I p 

i=l 1 < i < i 21 icl 


i=l 


( 1 . 2 . 1 ) 


Chung and Erdos [I 8, “] improve the bound (1.2.1). 
1.1 s»marl.es the reeult by On^g ^hcl 
hound is also reestablished by Whittle 

Lemma 1.1 f Sj - [l52^ 


n sf 

P ( U A . ) > — 

ViTl ~ S3_+^S2 ’ 

n 

provided that P( u A.) > 0 , where 

1=1 

n 

= 2 p . and Sr> = 2 -n 

1 < 1 < j in ir 


( 1 . 2 . 2 ) 


The lower bound (1.2.2) is exact if p. 
V 1 < 1 < j < n. 
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Lemma 1.2 gives the bound obtained by Dows on and 
Sank off Cwi], who improve the result of lemma 1.2, 13? emplo3ring 
linear programming technique. 


Le mma 1.2 10 21 

n 2s^ 2sg 

i (kTTi' “ E(En7 ’ 

where (k-1) is the largest integer less, than 
This bound is exact if n = 2. 


(1.2.3) 


2s,, n I 

and P( U A. ) > 0 ^ 
'^1 L =1 ^ ! 


Galambos Q 15^ shows that the lower bound (1.2.3) is 
the sharpest possible bound if the bound is to be e3q)ressed 
in terms of s^ Sg instead of 'P 

Gallot l?^] Kounias j2 29 2 ds'velop another kind of 
bound using matrix analysis. The following lemma summarizes 
their result. 


} 




Lemm a 1.3 [2 l*^!] » C U 
n m 

■ 2( u A.) > P . Q“.P, 

1=1 

where P = (P^jPg, • • • jT'ji) 0 Q” is an3^ 

of tna matrix Q = (P^-), satisfying = Q 

the f oil OY/ing elements 


(1.2.4) I 

! 

I 

generalised inverse i 

1 

. The matrix Q has I 


P.. = p. and P . - = P.. for all 1 < i i < n. 

XX id . d'i' ' ' . 

' ' ' 0 . 

In fact Gallot []ll?2 .obtains the above result for only 
nonsingular Q and Kounias (229 3 extends the result to singular 
cases. ' 
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Kounisis also shows that the bound (1,2.4) is exact if Q 
has rahfc one, since then all events coincide, or if it has rank 
tv;o i.e, there are essentially two events and one is a subset 
of the other i.e. 


El >P2 = or, Eg > Pi = Pi2 [.'29:1. 

It has also been saown by Kounias [^293^1 


n 

P( U A.) > Max 
i=l 



2 

i < 3 



1,3 e 


(1.2.5) 


where for 2 < r < n, {1,2, ...,n> having r elements. 


Sobel and Uppuluri give lower bounds of degree 

two for the cases when Pp = ^ ^ ^ 3 £ ri. 

Ihe lower bound obtained by them is 


n 

I’ ( U ) > 2iax 
i=l ” 2 < r < n 


r P. 


?n:il p 


12 


( 1 . 2 . 6 ) 


Phe most recent lov^r bound of degree two is obtained by 
Kounias and Blarin [l31)3 ♦ Phis bound provides an improvement 
over all other earlier bounds i.e. (1.2,1) _ (1.2.6), Kounias 
and Marin formulate the problem as linear programming problem. 
Since the present work is primarily based on their \;ork, we give 
details of the -work in following sections of this chapter. 

1,Z Lineav Vvogvamming Fo-mutat-ion 

In this section the Bonferroni lower bound problem is 
formulated as a linear programming problem, Phe following 
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notations will be used throughout. 

nCn+l)^^”^ consider the SL|±iI -dimensional yector space 

id . For notational convenience, we shall denote the 

components of any element x in this space a,s follows. 

n( n+1 ) 

file components of any X e R ^ are denoted by 

^i» ^ 1. Sind Xj^., 1 <, i < 3 < and are arranged as 

sho¥m below. 




- X, 


12 ^13 ^14 


^In ') 


--^Xgg - Xg^ 


X, 


2n 


9 


X. 


’n- In . 


n(n+l) 

Ihus by ^ G R we shall mean the vector 


z - (Xi,Xg,...,3^, Xj 2 ,...,Xj^,Xg 3 ,...,Xgj^ ... 

lo denote any vector In this space, it is only necessary to 

specify the element^., i < 1 < n, and , 1 < i <■ j < n. 

P e R 2 be the vector with the corresponding 


= P^ and xy 


•^13 • 


Similarly, f^^^wen J c H = 11,2,. ..,n, , we define 
the vector d(J) e H ^ “ as foliows. 

d{J) has components d^, d^., 1 < 1 < j < „ that, 
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arid 


= 1 if i e J , 

= 0 if i ^ J, 

d- • = -d. .d-, ¥l£i< 

Let denote the convex hull of the set {d(J) [ JC 1}, 

Let us consider the inequ8,lity 
n n 


P ( U A • ) > x_ + Z X,- 
i=l ^ ° i=l ^ 


^1 - 


1 < i < j < -n 


r p 

'ij iD 


(1.3.1) 

It is well lai own that, when x^ = 0 ohid x,- = x - • = 1,1 < i < j < n, 

^ 0 1 X3 ^ ^ 7 

the inequality ( 1,3 .1), which reduces to 

P( U A.) > 2 P. - 2 P.., (1.3.2) 

i=l ^ i=l ^ 1 < i < 3 < n 


holds for any arbitrary sequence of events 1 <_ i n. 

Ihe sharpest lower bound of the type (1,3.1) can be obtained by 
determining the variables x^» x . , x. 1 < i < i < n such that 
the right hand side of the inequality (1,3.1) tales on its 
maximum value. 

Let E be any sequence of e-vents A^^, A 
each e-vent is either the whole probability space or, null, Thus 
for such a sequence E, P- = 0 or,l and P. . = P- P.. Clearly 
d(.J ) defines the vector P associated v/itli the sequence E, when 
A^, i G J are full space and A^^, i ^ J are null. Such a 
sequence of events is called an elementary sequence. Obviously 
there are 2 elementary sequences., Renyi shov/s that 


2 ,...,A^ such that 
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any linear inequality of the type (1.3.1) holds true for an 
arbitrary sequence of events iff it holds for all tre 
elementary sequences, Kounias and Marin reestablish 

the result by using method of random variables, Galanbos 
Hailperin 19 3 also derive a similar result by using linear 
programming technique. Since P e 0^, Renyi’s result follows 
from the fact that if any linear inequality holds true for all 
the points d(J), then it also holds true for any arbitrary 
point of G^, 

Thus writing the inequality (1,3,1) for all the elementary 
sequences, we have 

Xq ^ 0, (1.3,3) 

X. + d(J)^’ X = + 2 X. ~ S X.. < 1, (1.3.4) 

° i e J i < 3 ” 

ij D e J 


f or all (f> J C H . 


Thus the sharpest bound of the type (1.3,1) can be obtained by 
Maximise 

n 

f y . . P , . C 1 ,5} 


c + I X. P- - 2. X. . p . , 

i= 1 1 <_ i < 3 < n ^ i ^3 


subject to the conditions (1,3,3) - (1.3.4). 
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Reform ulation of linear program ; 

Ihe above linear programming problem can be reformulated 
to a simpler form, witn the purpose to males the problem a little 
easier. Ror this reformulation the following results are useful. 

lemma 1.4 Csi!] 

If for any x^, z^, x^^, 1 <_ i < 3 < satisfying (1.3.3) 
(1.3.4), x^ is strictly negative, then the point x^ = 0, 

Xi X . . 

X- = , x-- = 1 < i < 3 < n also satisfies the se b * 

1 r_Xo 13 

of constraints (1.3,3) - (1.3.4). 

pro of : Since x„j x- , x- - , 1 < i < 3 < n is a feasible soluticn 
to the set of constraints (1.3.3) - (1.3.4) we have 

x^ + 2 *“ 2 X • • < 1 for any J { 1 , 2 , . . .,n} 

° i e ^ i < 3 ^ “ 

i,3 e 



{ _ 2 X. 

i G 


2 . 
i < 3 

1,3 S J 




==> 2 ±--2 x .-£ l . 

i G J . i 3 

ij3 e 

X. ^ X . . 

Hence the point x. = 0, x. = — , x- - = 1 ^ i < 3 < 21 

O' ' 1 X—Xq ■ 

satisfies all the constraimts (1.3.3) - (1.3.4), 



Thus for any solution of the system (1.3.5) - (1.3.4) 

such that Xq < 0, there corresponds another feasible solution 

(x , X., X..) with x^ = 0. The following lemma suows that 
' o^ X2 0 

feasible points of later type need to be considered for the 
optimization purposes. 


lemma 1.5 Csin 

There exists atleast one optiina.1 solution of tne linear 
programming problem (1,3,3) _ (1.3,5) such that x^ = 0 at this 
op t imal s olut ion . 

proof ; Assume that (x^, x'^, x*.) 1 <_ i < 3 < n, is the point 
at ?/hich (1,3,5) takes on its maximum value, subject to the 
restrictions (1.3.3) - (1.3,4) using the result of lemma 1,4, 
we construct a new feasible solution as 


Xo = 0, x^ = 


^i 

1-x! 


X 


ID 


. = 1 < i < j < n< 




The value of (1.3,5) at the point (0, x-, x-.) is given by 

^ ^D 


i ^ 

) i=l 1 <r 1 <’ -i ^ -n ID 3.3 


1 


z 


1 < i < 3 <, n 


( l-xj) ( 1-xp 


(1.3.6) 


n 


where z 


I P . X* - 


i=l 


11 


1 < i < 3 < n 


ID 'ID 


It is well knovm that, from (1.3.5) 
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==> X 


* * * 2 y r, 

- - ^0 S. ° 


^ *2 * * X 

2 - “0 - ^0 ^ ■<■ \ S . ^ 


=> (1-xJ) (x* + z) < z 




(1.3.7) 


The left hand side of (1.3.7) is the value of the objective 
function (1,3.5) at the point (x^j x*, ^c*^). Hence^ if 


__ z_ 

1-x: 


10- 

(x*, xf, xf^) is the maximum point then, + z = , 

Hence (0, x- , x-^) is also the optimal point. 

id- X J 

This completes the proof . 

Thus using lemma 1.4 and 1.5, the variable x^ can be 
dropped from the linear programming problem. Hence the resulting 
linear programming is : 


n(n+l) 


I*P(I) J 


To determine X e R 


n 


, which 


llaximizes 2 - - 2 P.-x-^ (1.3,8) 

i=l 1 < i < 3 < n 

# " 


subject to 2 x^- — 2 x - ^ < 1 

i e J ‘ ^ i < 3 ^ ” 

ijj e J 


(1.3.9) 


for all nonempty subsets J of If* 

We will no?/ onwards be concerned with the problem liP(H). 
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1.4 Mot-ivcction 

Problem IP(l'T) is' a linear program of ^ unrestricted 

■variables subject to 2^-1 constraints. As n increcises, the 

number of constraints becomes unsurmouiita.ble , end nence none of 

the the conventional linear programming algor itnms is capable 

of effectively handling such large size problems. However some 

of its special features are encoura.ging . It is observed that 

for fixed n, the set of constraints (1.3.9) remains independent 

of the sequence of events of their respective probabilities 

P :| j P^ ■ and so does the convex polyhedron defined by the sex of 

constraints (1.3.9). Moreover, for different values of n the 

coefficient matrix of the set of constraints being similarly 

structured, the geometry of the con'vex polyhedron is similar 

n(n-fl) 

and hence the set of extreme points of the polyhedron in R 

may follow a particular structure. This important feature of 

LP(IT) is due to the special structure of the constraints and of 

the right hand side constants. Hence, notv/ithstanding the size 

of the problem, it may be possible to develop a special purpose 

algorithm vi^hich can effectively exploit the structure of IP(N), 

■» 

Motivated by this observation, an attempt is made to explore the 
structure of IP(R) and to develop an approximate algorithm based 
on its special structure. 

2.5 Basic Approach 

The approach followed in the present work is mainly 
influenced by the work of Kounias and Marin ^31^]. This section 



13 


gives a brief outline of their approach. 

Let I’(II) denote the convex polyhedron defined by the 
set of inequalities (1.3,9). denotes the subset of N, having 
h elements ¥7here N = { 1, 2, . . . ,n > 

Ite fiii it i_cm _ 1. 1 

X e R ^ is said to be an extreme point of i'(h) 

if it satisfies linearly independent constraints of 

(1,3.9) as equalities. 

Obviously, the convex polyhedron l'(lT) is unbounded. But 
since the maxirum value of (1,3.8) is bounded ao 1, the maximum 
value will be finite and hence it is attained on atleast one 
extreme point of F(N). In a rather unconventional way, ICounias 
and Marin define a point X G F(K) to be a vertex point of P(N) 
if X satisfies constraints as equalities and these 

constraints ne ed not be linearly independent. Since any extreme 
point is also a vertex point, it can be directly concluded that 
atleast at one of the vertex points the objective function talces 
on its maximum value. Ihus maximizing the objective function 
(1.3.8) over the set of all vertex points of I'(ii) V7e geb the 
optimal value of liP(N). 

The approach used by Kounias and Marin is to generate 8 
classes of vertex points, V 2 ,...,Yg, where the optimization 

can be carried out. Though, in the process, all the vertex 
pointsfor n <_-4 are obtained, these 8 classes do not cover all 






14 


the vertex points of ) for larger values of n. However the 

8 ' ■*' 

ii^aximum value of (1.3.8) over \J V- gives a lower hound letter 

i=l ^ 

than any of the previous bounds discussed in section 1.2. 

Ihese 8 classes of vertex points are as follovTS : 


Yi : Xi = 


^ X* , i e J^, 


* 


0 


i ^ J, 




■13 ’ ^’3 e J^, 
i < Dj 

0 5 otuervy'ise, 


3 < r < n, Vi/here x^, Xj_^ are the coordinates of any vertex 


point of ^(J^) 




(1.5.1) 


Vg : Xf = 


1 S 
) ^ 


5 i e 




) ijD 6 


j kTk-iT 

(^ 0 , otnerwise, 

for all 2 < k < r_l and k = 2, if n = 2; (1.5.2) 


• ^ 1 

0 , 1 ^ ^ 


2 

E . i e J 


¥3 : 


n-l» 


"^ii = 


j- ® ^ 1 * 

S[ 

where J^f) = <|, , for all 2 < k < n-2 


\ k^kTiy > iyj e 

2 


(1.5.3) 


Y 4 : X 4 




r 


V 


X - - = 

3-a 


£ j i e j 


i» 


k^k-lj » ij3 S 

1 : ■ . • . 
kXKZTJ J je 


V 


^1^ ^n-1 = ^ y 3 1 ^ 


9 


(1.5.4) 
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Yg : Xj_ _ 


1 f j ^ 

I 


’ ® ^n-1’ 




^iD 


1 


v_ 




" ETk'^T’^ e j^sD e 

Jl n =<p , for all 2 < k < n_3; (1.5.5) 

2 


r 

t 2 


E 


e J 


n-2’ 


^6 = = / 


^ij “ / 1 


"ICT ’ 1 e Jg, 


\ ^ ^n_2» 

kpiTj^ i e Jg, j e 
0, 1} D e J2’ 


f or Jg n p = 4) j fo^ 2 < k < n-Sj 


(1.5.6) 


^ i > ^ ® Jn-8> 


Vn • — 




k(k-fy 


5 i, D e J 


n-2’ 


k 


, i e J 


2’ 


‘iD / 1 . 

'\ klklTJ ’ ^ e Jg, 3 8 


0 1,3 e J, 


V 


Jp 0 2 ~ 2 ^ 5, 1^-2^ 


(1.5.7) 


V 


g : = 2bj _ -b^ , Xy = Sbj^b., 1 < i < j < n, 


(1.5.8) 


wbere Z b . = 1 for atleast SlSiil subjects J of 

i e J 


Ihe following results are due to Kounias and Mar in []] 31 "I 

Let be the maximum value of (1.3.8) over X 6 for 
1 < i < 8. 

Theorem 1.1 

n 

PC U ) > Max L. 

L=1 “ 1 < i < 8 ^ 


(1.5.9) 
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Theo rem 1.2 

Max L- is sharper lower hound than the hounds 

1 < i < 8 ^ 

( 1 . 2 . 1 ) ~ ( 1 . 2 . 6 ). 

Thus Kounias and Marin give the best knovn'i lower hound hy 
maximizing (1.3.8) over the 8 classes of vertex p oints . But it 
is not the best possible hound always. Moreover for the 
maximization purposes it suffices to concentrate only on the 
set of extreme points, instead of the set of vertex points 
which is larger than the fomer set. Keeping this in mind it 
can he said tha.t a successful characterization of the set of all 
extreme points together with a tractable algorithm to optimize 
the objective function over the set, would give a better v/ay of 
computing the best lower Bonferroni bound than any of the 
conventional approaches. A partial characterization of the set 
of extreme points is given in this thesis. 

Follov/ing the same approach of Kounias and Marin, a 
broader class of feasible points is def ined. The maximum value 
of the objective function over this class improves the bound 
obtained by Kounias and Marin. Since the trivial way of getting 
the best bound is to solve the linear program by simplex method, 
though at a high cost, emphasis is given on nuaierical computation 
of best lower bound at a low cost* For several special cases 
some simple v/ays of ccmputing the exact lov;er bound are discussed. 
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1.6 Thesis Summary 

Besides tiie current chapter, the thesis contains 4 
more chapters. 

In Chapter II we introduce a nev; class of feasible points 
V* of P(N) . It is shov\n that the maxhiium value of (1.3,8) over 
this class provides the present best lovjer bound. Since, no 
efficient technique to compute such improved bounds could be 
developed, give tv7o different computationally trantable 
algorithms to compute two different types of lov^er bounds. It 
is also shown that one of these two bounds is alv;ays better 
than previously loiovn bounds, numerical computation of improved 
bounds involves minimization of a quadratic convex fimction of 
0-1 variables. A branch and bound technique is proposed to 
solve such problems . Besides finding lo^^er bound, attempt is 
made in this chapter to solve linear program LP(E) through its 
dual by the use of column generation technique. It is shown 
that subproblems arising from column generation teclnique are 
equivalent to minimization of a quadratic Pseudo-Boolean function. 
Computational aspects of such problems are nor explored to its 
full extent. 

In Chapter III, we deal with the structure of tne extreme 
points of the convex polyhedron i'(N), Iv^o different kinds of 
extreme points are characterised. Ihe first kind is characterised 
in terms of the coordinates of the points, Por the second kind 
of extreme points, though the explicit characterisation of the 
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coordinates is- not given, it is sho\m that such points can he 
obtained from a type of feasible points whose coordinates are 
known. Besides this, a partial characterisation of all the 
extreme points of E'(N) is given. Every attempt to generate all 
the extreme points of E'(N) proved abortive, however, some 
developments in this direction are reported aaicl a conjecture is 
stated leading to generation of all the extremo points. The 
adjacency structure of special kind of extreme points is also 
discussed in this chapter. 

In Chapter I¥, some particular cases are studied for 
which the best Bcnferroni bound can be computed without much 
effort. 

In Chapter V, we explore the possibility of generalizing 
the methods studied so far, for the cases of Bonferroni bounds 
of higher degree. We show that the class 1 * can be generalised 
for such cases. The structure of some types of extreme points 
of the corresponding polyhedron is also given. 




Chapter II 


BBS I LOM BOTm i> 

2.2 Introduction 

Ihe problem of best Bonferroni lower bound is formulated 
.in section 1.3 as the linear program : 

lo maximize 

z= £ (2,1.1) 

i=l ^ ^ 1 < i < 3 < n 13 

IP(H) : 

subject to £ X- — £ x. - < 1 (2.1.2) 

ieJ i<D 

e J 

for all none irpty subsets = {l,2,...,n>. 

Considering the difficulties encountered in solving IP(B), 
attempt is made to maximize the linear function (2.1.1) o^r a 
smaller set of feasible points. The maximum value of (2,1.1) 
over any subset of tiB feasible region provides a lower bound 
which may not necessarily be the sharpest bound, A desirable 
characterization of such a subset of P(l) is that it should be 
small enough to be easily handled for maximization purposes and 
large enough to give a considerably stringent bound. In this 
chapter a class of feasible points, denoted by Y*, is defined 
and it is shown that tie maximum value of (2.1,1) over this 
class is sharpest among the knovoi boun(|s (1.2.1) - (1,2,6) and 
(1,5,9). The maximization problem over the class leads to 
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maxiniizitig a segmenoe of unconstrained qaadratic functicais of 
n integer variables* It is to be noted that, though the liP(N) 

is of variables, a reduction of tbe number of variables 

is obtained at the expense of the linearity of the object tire 
function* Similarly, it would be seen that by imposing integer 
restrictions on the variables, the set of constraints (2*1.2) 
are automatically taken care of. Though no exact method to 
solve such quadratic- integer optimization problem is discussed, 
an approximate algorithm to get a lower bound is developed by 
defining a type of local maxima* Such approximate method leads 
to maximzing a Pseudo-Boolean quadratic concaire function. A 
branch and bound technique is develflpe^ to aolve such a problem. 
We report our computational result in section 2.8. In section 
2,9 we show that in order to get the best loiown lower bound one 
needs to maximize a concave Pseudo-Boolean quadratic function 
of 2n variables. Such problems can be solved by the proposed 
branch and bound method. In the last section we atteE 5 )t to 
solve the dual of IP(If) by column generation technique and show 
that the subproblems are equivalent to minimizing Pseudo-Boolean 
quadratic functions. 

2,2 The Clc&B V* 


In this section a class of feasible points is defined and 

some properties of this class are studied. 

n(njt-l) 

Let the class Y*C. R ^ be defined as follows t 


Por any X Q Y* 
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(2h+l) - a| 


1 < i < n 




fi 


h( h+l) 


1 < i < 3 < n 


( 2 . 2 . 1 ) 


where a^, 1 <. i < ii> snd. h ^ -1,0 are arbitrary integers. 

lor any given (n-i-1)- dimensional integer vector (a, h) = 
(aijag, .. ,,anili)j using (2.2.1) a point Z G Y* can be obtained. 
There can be several such integer vectors (a,h) which define 
the same point X 6 V*. lor instance, 

let a^ = h, ag = h+1, a^ = 0, i > 2, h arbitrary. 

Then the corresponding X is 


= 1, Xg = 1, x^g =2, XjL = = 0, for 2 < i < j < n, 

( 2 . 2 . 2 ) 

Clearly the point (2.2.2) is generated for any arbitrary h and 
hence there are inf in it ely many (a, h) v/hich define the same 
point X of. Y* • In tteorem 2.1, we shov/ that we can restrict h 
to be a positive integer and still get all the points of V*. 


Theorem 2.1 


Any X G ?* obtained by (a, h) can also be obtained by 
(-a,-h-l) using the relation (2,2. 1)^ 

Proof : Ircm (2.2.1), we have, 

' ' '2 ■ 

’'i “ “ RbO)' • 1 < i < 3 < n. 


The point generated by the vector (-a,-h-l) is given by 
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{2(-.h-.l) + (2in-l) aj_ - a| 

(-h-1) (-h) ~ h(l3+l) “ 

2(-a^) (-a^) 2a^a-j 

This con^letes the proof* 

Thus, in order to define all the points of Y* it suffices 
to consider the (n+l)-.dimensional integer vectors (a, h) with h 
positive* Although there is no one— one correspondence between 
such set of integers and the points of Y*, only for the sake of 
simplicity we shall notationally use (a, h) G Y*, h > 1, as the 
point X Q Y* obtained frcm (a, h) by (2.2.1)* 

n(n+l) 

^ p 

Let Y-^ denote the class of points X G E. obtained 

from (a, h), h >_ 1, by the fomula (2*2*1) where the integer 
restriction on a is relaxed* Similarly, we define Yg as tte 
subset * of Y* where in relation (2.2.1) a^’s are restricted to 
tl^ values —1,0,1 and 2, and h > 1, integer. For a fixed h, 
define a subset Y^(h) of Y^ as; 

Y*( h) = {(a,h) I a is an integer vector} * 

The subsets Yj^(h), Vg(h) can ,|)e similarly defined. Let 

L^, L*, L^, L^(h), L*(h), I^(h) be the maximum values of (2.i.l) 

over the classes Y^, Y*, Yg, Y^(h), Y*(h), Yg(h), respectively. 

The following theorem shows that Y* is a subset of the 
feasible region F(l!l)* 
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■Theorem 2.2 


V*OI'(K) 

Proof ; Let us consider the constraint (2.1.2) corresp ending 
to any arbitrary subset JON, 

Por any X = (a, h) e 1*, substituting the values of x-,x. . 
from (2,2.1), we have. 


2 X. - 2 X. . 

i G J i < D ^ 


= 2 
i e J 


i,3 e J 

( 2h+l) a^ - 
h(lH-l) 


2a. a. 

2 — ■ J. 

i < 3 h(h+l) 

i,3 e J 


h( hf 1) 


(2hfl) 2 a. - ( 2 a.)' 

iGJ ieJ 


h( ht-l) 


Fh - 2 a- "IF 2 a. - h - ll + 1 

1- iGJ -Jl-iG J -J 


( 2 * 2 . 2 ) 


Por integers a^ and h ^ -1,0, tte following inequality holds. 


i — /h- Z a. } {*£ a..-h_l)+l<l (2.2.S) 

ieJ^ leJ^ 


hence, from (2,2,2) v?e have. 


2 X- 
iGJ 


i < D 
i,D G J 


< 1 
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Since the inequality (2«2,3) holds for any JC N, hence 
the point z = (a, h) e V* satisfies all the constraints (2.1,2). 

■“ Hence Y*CZ F(N) and this completes the proof. 

In the theorem 2.3, it is shown that the class Y* is a 
more general class than any of defined in [][3l2] . Let 

convex hull of V*. 


fheorem 2.5 

for 2 < i < 7. 

proof j We prove the theorem in parts. 

(i) In this part it is shown that ^ 2^1 Y*, Y^C Y*. Putting 
aj_ =s 1, i 6- J and a^ = 0 i ^ J in (2.2,1) for any J ON, we get 


? ? 

X. = ^ i G J, X.. = 

^ in-1 13 h(h+l) 


X,- = X. . =0 otherwise , 
1 ID 


i < D» i »3 e J, 


Prom (1.5.2) clearly f or h = k-l, 2 < k 1, 1*^1 ~ ^ 

if n = 2, the above point is in Yg. Substituting a^^ = 1, 
i e -a^ = -1, i e in (2,2.1) we get 


t 


■'I' 


X. 


1 hfl ’ 13 h(hfl) 


i i i, i»3 e J. 


IL^l 


x^ = - 


h 


^3 h(hfl) 


ie Jj, 3 e 


Prom (1,5,3) clearly for h ss k-l, 3 £ h < n-1, ^ > 

the above point is in Y^. 

(ii) In this part, it is shown that points in Y^, ^ < i < 7 
can be eapressed as the ccaivex combinations of points in Y*, 
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let us consider tte pair of points (a’,b) and (a", h) where, 
a5_ = 1, i e 

a[' = 1, i e Jn- 1 » ^-i = ^ ® '^n -1 = <*' * 

Then the point X 6 V^, from (1«5,4), is gi^ven by 

X = |x* + I X" where X* = Ca’,h) e V* and 

X" = (a”,h) e V*. 


2 _j 2 ^ ^ ^ „ -r 


n 


To see this, from (2 ♦2*1) we have, 

~ ^ * ^iD h(krf)' 

^i = F » ^ij = ic(k-iy ^ i»D e 

♦ 

and x^' = 0 = x”^ for i G j 6 


Thus Xi = | x’ + I xj = J 


r 2 - j 

k » ^ ® ‘^n^l 


I , i e 




r 2 


\ k(]&-l) ’ 


1 i 1 « 

^iD =2 ^i3 + 2 ^i3 = / 

^ 1 


i < j» i,3 e J. 


n-l» 


k(k—l) 


, i e J 3 _, 3 e 




Thus from (l.S.d), for 3 ^ k < n-2, the point X is in 

Similarly, any point X 6 Tg can be e 35 )ressed by the points 
X*j X” 6 Y* defined as follows; 

X* = (a’,h) G Y% such that a? = 1, i G = 0» ^ 

X" = (a”, h) G Y*, such that a” = 1, i G si” 5 = —1, i G 
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where ^ correspond to subsets defined in (1,5.5). 

It is merely a routine to check that for 2 <_ k < n— 5 any Z Q 
can be written as, 

X = I X* + I X". 

Similarly the following pair of points give rise to the point 
in Vg as their convex combination, For any ^ S ¥g, define 

x' = (a’,h) G Y* such that 

a!^ = 1, i e = 0, i G , a^L = 1» ^ 

and x” = (a”, h) G Y* such that 

a? = 1, 1 G a!* = 1, 1 G aj = 0, i G 

where J 2 = JiU and Jj^2 defined in (1,5,6), 

Thus for h = k-l, 2 ^ k < n-2, any X G Vg can be written as 

X I Tr ^ , 1 TT ^ ^ 

= ^ ji. + ^ X , 

Similarly for the case of X G Yrj>, i,e, from (1,5,7) 

^i = eIt > ^iD = fc(k;i) * ^ ® 


H = E+T » ^ ® ‘^2» ^ii “ 


i3 = k(^ * ^ ® 3 e <F^2»^i3 = 


we define the pair of points X* = (a*,h) G Y*, X® = (a”,h) G Y 
such that. 


= 1, 1 e Jji_2i - 1, 1 e aj _ 0, l 6 
aj = 1, i e J^g, aj = 0, i e Ji, a" = 1, i e 
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Then, clearly 1 = ^ Z". 

Thus part (i) and part (ii) complete the proof* 

Oorollary 2*1 : 

G 2 < i < 7. 

Proof t It is to he noted that in the proof of theorem 2.3, 
in order to show that "7^^ O 0 V* for 2 < i < 7, the -values 
required for any are o, 1, -1 only. Thus for any X G 
2 < i < 7 there exist X\ Z“ G Tg such that X G G[]]^ X* ,X**]]| • 

2* 3 Best tower b o un d 

Clearly V*^ P(N) 

==> Ii* < Max of (2.1.1) subject to (2.1.2) 
n 

==> Ii* < P( u A.). 

i=l ^ 

* n 

Thus 1 is a lower bound of P( [j A^) • In this section, it is 

i=l ^ 

shown that i* iiaproves the bounds obtained by Kounias and Marin 
^ 31 and hence from tteorem 1.2, 1* is tte best known bound. 

Theorem 2.4 

> Max Iij * 

“ 1 < i < 8 ^ 

Proof : In view of theorem 2.3 we need only to pro-ve that 

L* > Max iH> %)• 

The prooJ^ is in two parts, the first part shows L* > Lg and the 
second part, L* ^ L^. 
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( i) Substituting tte value of Xj_, Xj_^ frcm (1*5,8) tte 
objective function is, 




Z 

1 < i < 3 < n 


2b, 


’=3 


3-D 


= 2 P^b — b^Qb where b = (b^,*..,b^) 6 (2,3,2) 

Ihe function (2*3,2) takes its maximum at the point which 
satisfies 

Qb = P 

or, b = Q"*P for a generalised inverse Q” of (2*3*3) 


The maximum value of (2,3.2) at (2.3.3) is given by 


2P^Q“*P - bVQ“Q’b = 2P'^Q*1> ~ PQ“P = PQl? (2 *3.4) 

Since Yq is defined by b^, 1 ^ i < n from (1,5,8) where b^^s 
are subjected to restriction that at least for constraint 

Z b^ =: 1, hence clearly 

# 

Ig < P^Q~2 

Since any irrational number can be approximated to a rational 
number to any degree, it may be assumed that b^*s are rational 
numbers. Thus, 


where p 
h = ICM 


bj^ = ~ for all i, 
^i 


. , q. are relatively priioe numbers and q, > 0. 

U: . wim JL. 

P-f 


b. 

i 


^i 


(2*3.5) 

let 

(2.3.6) 
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v^here is integer, for all 1 < i < n. 

Let us consider the point (a, h) 6 V*, where a^, 1 < i < n and h 
are defined in (2,3.6) • 

Substituting the value of (a, h) G T* in (2,1,1), we have, 

2 

n (2h+l) a- - a- 2 a. a. 

V - - - - P , 2 , - - A _iL p ^ 

i_l h(ln.l) i l<i< 3 <n 


= {(2hfl) - a\a> (2.3.7) 

Substituting a = b from (2,3,6) in (2.3,7) 

= {(2hfl) P^b - h b®Qb> 


= {(2hn-l) P^Q"P - hP®Q~P> 
= P^Q“P, 


(ii) In this part we show that L* ^ L^, 

n(n 4 -l) r( r-f 1) 

Like X e R ^ , let X^, a point in E ^ , be defined as 


= ( 

^ ^2 ^ ^ ^ ^12 ^ ^ ^ lx ) X(0,r) denotes the 

Si£+il 

point in E obtained frcm as follows j 

x^(0,r) = Xp 2c^^(0,r) for i <3, i,3 e N. 

Xj_(0,r) = x^.(0,r) = 0 for i 0 

Si£p:l 

Let V* be the set of points of the kjthd V* in R l.e, 

^ defined by r+1 integers a^, ie J^, h > 1 frcra. 
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the relation (2.2.1). Similarly T»e define as the vertex 

p oints of the kind in Clearly, 


e Y* ==> X(0,r) G Y*, 

and the value of the objective function at X^, an<l-X(0,r) are same. 

Similarly, X^ G ==> X(0,r) 6 ¥ i = 1,2,. ..,8. 

Ihe class Y^, from (1.5.1), is the class of points of the type 
X(0,r), where X^ is a vertex point of P(Jj,). Since ttere are 
only 8 classes of vertex points known so far, we shall use the 
following definition of for all practical purposes 

Y^ = {X(0,r) I X^ 6 ij (V^)^ for all 1 < r < n >. 

2 


Hence 


s= Max Max \ 2 ^ 

^ 1 < r < n~l,J^ - 8 I i G ^ ^ i < D 

- “ X^ G U(^i)r 


5=2 


i>3 e J, 


V. 


Max 


Max 


2 2 

1 i i \t Ij i;| 


1 < r < n-l,J^ ^ g (V*)„ l ^ ® '^r 


I 


ij;i e d, 


n 


Max 


Max 


^ ^iX-f(0,r)— 2 P- ■x..(0,r) 


1 < r < n-l,J^ X(0,r) 6 Y* ] i=l ^ ^ ' l<i<j<n f3"i3’ 


< 


Thus 1^ < L*. 


This completes the proofi 
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Corollary 2«2 

2 < i < 7. 

Proof : Proof directly follo?;s from coroll^iry 2*1* 

==> I* > Lj_, 2 _< i < 7. 

2,4 Numepiaat amputation of tcwev hound 

In this section computational aspect of the bound 1*“ is 
discussed. It is shovm that the problem of computing 1* leads 
to hard combinatorial optimization problem. 


Substituting the value of X-- from ( 2 . 2 , 1 ) in (2,1,1) 

i X J 


T/ve have 


Max zl = 


(2hfl) - a^ 


Qa 


h( h+ 1 ) 


(2.4.1) 


PI : 


subject to a G ^5 h > 1 integer, Y/here Z is the set of 
integers . 

Lemma 2,1 

For a fixed h > 0, zl is a concave function of a for all 
a G . 

Proof : For fixed h, 

L (2^^) ^ “ a^Qa]] 

is obviously a quadratic function of a 6 R . Since Q is a 
covariance matrix d 29 3 » is positive semi-definite. So zl 
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is a concave function of a e pA Hence tHo pnoOfen of oo^utln, 
(a; reaucee to maxtalzlns a concave cuadnatlc function of n 
integer variables Ihe apparent sniplieity of the problem Is 

deceptive, ge fall to develop anv effc„+- 

P J ifective al^.oritlim to solve 

such problems. Por some unlmovvn ree=on. such ai 

xud^ooiib sucn 'oroblerns are 

constantly ignored even In literature. Some rerevonit references 

“ V C30, L36D 

psj, L-: - nab:,. Basically tno types of probj 

a^ been dealt mth so far. Pirst, toe maximisation of any 
omstraioied quadratic function of integer variables 
C 9 J > E 21 2] ~ 28 “] , r 42 “1 _ r 44 -I 

'i- -J • Second, the maximizatio] 

of concave quadratic function of int-i^er .r r • 

. J-otc.ger Variables subject to 

ine^^ constraints and non-negative restrictions on the variables 

I- J - CSJ, LSZZI, ISBJ, [;36.-I, E583. It is found that 

none of these algorithms are suitable -i-n +• i i 

uixaole to talce due advantages of 

the structure of the problem -o 

„ , ^ ^ ^ consideration. She algorithms 

0 the first kind are incapable of makine: n- +vo 

^6 use oj. the concavity 

of the function and so thev are nn+- wa 

y are not recommended for the present 

problem. Ihe second type of algorithms are he-vilv de c a * 

°-oe necvvily dependent on 

the set of constraints, that is on th« v , 

J IS, on the boundedness of the 

feasible region. Hence these alvorithmo a. 

algorithms are not efficient in 

solving any unconstrained maximization problem. 

Thus it is difficult to develop an +• 

^ V lop an exfective algorithm to 

f ind 1 (h). Moreover, l^(h) which dene-n^t, 

V J wiicn aepends on n, camiot be easily 

expressed as a simple algebraic function of h. Unless the 
property of l*(h) as a function of h la even nn the 
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presence of effective algorithm to maximize unconstrained, 
concave quadratic function of integer variables , it may he 
difficult to determine 

L* = Max l*(h), 
h > 1 ■ 

h integer 

In view of these difficulties an approximate method to 
solve PI is given here. The approximation of L* in the proposed 
algorithm is made with respect to the following conditions, 

(1) Instead of searching for a global maxima of concave 
quadratic function with integer variables, a kind of local 
maximizing point is searched. This local optima is obtained 
by rounding off the global maximizing points when the variables 
are not restricted to be integers, For the problem under 
consideration such approximation seems to be good, Por all 

the problems encountered during computation such local maximizing 
points were also observed to be global. 

(2) In order to compute L* = Max I)'**(h), h is restricted 

h > 1 

to be less than n. This assure t ion is based on the coij>utational 
experience. Thus there are only n quadratic maximization 
problems to be solved, , 

With the help of these assumptions an algorithm is 
developed which finds an approximate value of 1* in a finite 
number of steps. 

The locally maximizing point is defined as follows; 
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Definition 2,1 ( Locally maximiziii^ points ) 

let a* lie the real maximum point obtained by relaxing 
the integer restriction* let the components a* of a*, i G J, 
be fractions and for i ^ J be integers, v/here |j j = m where 
^ 1 ^ ^ locally maximizing point is the best point among 

the set of 2^ integers obtained from a*, by fixing the 
conponent a*, i ^ J and by rounding-up or, -down the components 
S.*, i e J • 


Yife concentrate on finding such locally maximizing points 
of Pi, This problem, as we show below, boils down to a 0-1 
optimization problem. 


Por a fixed h, a* can be determined by method of 
differentiation a^, 

a* = 2^ (2.4.2) 

Thus substituting the yalue a* from (2.4,2) we haye, 


l*(h) = 


(Eh+l) a*P - a*®Qa* 
h(li+I) 


2a*®Qa* - a*®Qa* 


h(h+l) 



(2.4,3) 


Without loss of generality, we can assume that all the components 
of a* are fractional and let 6 j|_ be the fractional part of a"^. 

So we have , 
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~ n ^iZl (2 *4 ,4 ) 

where |a*3^ largest integer less than a^* let 6 le the 
vector (6^, S2 ****»^n^ S • For such a a^ and hence for a 
given 6 we define a set of vectors 

W = t a* - 6 + y I y e Bg > (2.4.5) 

H 

where Bg denotes the set of n-dimens ional vector with coinponents 
0 or, 1, 

Row the locally maxinizing point is obtained by solving 

Max (2h+l) c^ - c^Qc (2,4.6) 

c e w 

Substituting the value of c for (2,4,5) in (2,4,6) 

we have 

(2h+l) (a*-,6+y)^ - (a^-6+y)'^Q(a^-6+y) 


= ^(2h+l) (a^^~6^+y^) - a^^Qa* - 6^06 - y^Qy 


+ 2a^Q6 + 26 Qy - 2a^ Qy 


= (2h+l)a*^-a'*^Qa*-6^Q6-y^Qy+26^Qy 

~ (2h+l)6^ + (2h+l} P®y + 2a*®CB-2a*^Qy J 
(using relation (2,4,3)) 

= (2h+l)a*^ - a*®Qa^ - 6®Q& - y^Qy + 26^Qy 


= M - y^Qy + 26 ^Qy 


(2.4.7) 


% 
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where M = (2h+l) - a*®Qa* - 6®Q6 . 

Thus the problem (2.4,7) is equivalent to 

P2 : Min y'^Qy - 26^Qy (2.4.8) 

y e b| 

Let Zj(h) denote the optimal value of P2. 

2,5 Algorithm 1 : to solve P 1 

In this section ?,e describe the proposed algorithm to 
get an approximate solution to the problem pi. 

Alg or it hm 1 : To solve PI, 

Step one : Initialisation! h = 1, u = n 

step two : . Pind a'^ = q“p, I^(h) = P'^Q“P. 

Step three ; Solve the problem P 2 for the a* £ derived in 
step two. The voAue Zj(h) is computed. 

Step four : The current Bonferroni lovrer bound is computed as 

L = (2h+l) a*^-a‘*^Qa*-6^Q6 - Zj(h). 

Step f ive ; Compute the integer k such that 

i 

St ep six .; Update u = Min (fc,u) 

h = h+1. 

S tep seven ; If h > u, stop, otherwise, go to step t\;o, ' 

The efficiency and finiteness of the algorithm depend on 
those of solving the problem P2« But for the step three, the 
algorithm is simple and has at most n steps. 

* ■ . ■■ 



37 


Thus the whole attempt to get an improved lower bound 
boils dovai to solving the problem P2. In the follovjing section 
we shall discuss the techniq.ue to solve P2. 

2.6 Branch and bound algorithm: for general 0-1 problem 

The problem P2 can be viewed as minimization of Pseudo- 
Boolean function which is convex and quadratic. The importance 
of O-l program has led to the elaboration of numerous techniques 
for their solutions'. Techniques which deal v/ith optimization of 
quadratic convex polynomial are mainly branch and bound or, 

Boolean method and are discussed 21 '] - L 28 3] j " 

C I] > E^^El* these approaches are iterative, 

each iteration consisting of two phases, one of analysis of the 
given problem such as fixing the variables, finding good bounds 
for optima, and the other of synthesis of the new subproblem by 
branching or, back tracking. Boolean method is found to be 
more useful for analysis rather than for synthesis. Boolean 
analysis is not helpful for the problem P 2 . T he p r op ose d 
algorithm uses analysis similar to that used in works by Hammer 
[]21[3 “ except that it does not malce use of Boolean 

analysis. On the other hand, it follows the usual branch and 
bound technique suitably modified for the structuie of the 
problem P2, 

We give here, a brisf outline of the branch and bound method 
to solve a general o-l optimization problem and develop the same 
with reference to P2 in section 2.7. Notation for branch and 
bound is followed as given in Murty E E • 
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Branch and bound method provides a methoc-olory to search 
f or an opt imum feas ible solution by enuunrating partially a set 
of solutions and implicitly eliminating large groups of potential 
solutions to the problem without ejqplicitly evaluating them. In 
the course of applying the branch and bound approach, the overall 
set of feasible solutions is partitioned into sma.ller subsets. 

At each stage one subset is chosen and an effort is made to 
find the best solution in it , If the best solution is found, 
then the subset is said to be fathomed.. If it is not fathomed, 
then the subset is again partitioned into tv/o subsets and the 
same process is repeated, Ihe problems associated with the 
unfathomed subsets are called can dida te pr obl ems. 

With reference to a 0-1 optimization problem, the candidate 
problems are obtained by selecting a subset of variables and 
fixing them at 0 , or, 1, Any variable that is fixed at 1 is called 
l~variable and any variable that is fixed at o is called 
0— variable , The 0-variables at 0 and l—variables at 1, constitute 
^ p8-htial s olution . Each candidate problem generated in the 
algorithm corresponds to a partial solution, Yariables that are 
not included in the partial solution are called free variables 
in the candidate problem. Given a partial solution, a completion 
of it is obtained by giving values 0 or, 1 to eacn of the free 
variables, 

We ill first discuss the analysis to be performed cm a 
typical candidate problem. Consider the candidate problem in 
which Uq, U^, are the sets of subscripts of the 0-, 1- and 
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free -variables respectively. iPor the given candidate problem, 
we shall determine tests, 'ISSll , '1BS12, _1EST 5 so theit lES'll 
gives the fathoming criterion, 2ES12 helps to augment the 
candidate problem to obtain a new candidate problem, and 'J3SS13 
partitions the set of potential solutions to t',,’o different subsets 
and defines two different candidate problems. If the current 
candidate problem is fathomed, a nev/ candidate problem is picked 
up from the list of unfathomed candidate problems. Here, the 
candidate problem which is the latest entry in the list is taken 
as the current candidate problem. If the current candidate 
problem is augmented to a nev; candidate problem then the new 
candidate problem is taken as the current problem. If two 
candidate problems are generated from the current candidate 
problem then one of these is taken as the current problem and 
other is entered .into the list of unfathomed problems. 

At any stage, i ncumbent is the feasible solution that has 
the least objective value among the feasible .solutioi.'S obtained 
from so far fathomed candidate problems. Let LL be the lov/er 
bound of the value of the objective function wirh respect to 
the current candidate problem. If LB is greater than the value 
of the present incumbent then the candidate problem is deleted 
from the list and this operation is called pruning, lor pruning 
purpose, the lower bounding strategy ?i?hich we adopt here, is as 
f ollowsi all the integer restrictions of the free variables are 
relaxed and the minimum value of the objective function for the 
relaxed problem gives the lower bound, LB, 
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2,7 AtqoT^tlm 2: to soli^e P2 


• We now proceed to discuss the main features of the present 
algorithm with reference to P2, 


She general form of P2 is 
^ n ^ 


1 - ■ y ■ y • 


Minimize -2c.y.+ ZP.v'^ + 2 L 

i=l i=l 1< i< 3 <n-^^ 

(2.7.1) 

for y^ = 0, 1 for i e where c^ = 2(&Q)^, defined in (2.4.8). 
The typical candidate problem can be written as? 


2 2 |c.-2 2 p. 

i e ^ i e 1 D e % 


Wi + 2 p 

i G ^ 


+ . ^ 2’i yp + 2 2 y. y • + 2 2 p . . 

leu^ i <3 13 ID i<j 

ijD e U£ i, j e b\ 

(2.7.2) 

= H(U^) 2 c.(U.) y - + 2 p.yf + 2 2 P-. y- y 

-oTT -J- jr^T-x i- i --- 


1 e i e 


i < 3 '^ 3 - 3 

1 } D S 

(2.7.3) 


whe re 


H(U^) = ~ 2 c . + 2 P - + 2 2 

^ O TT -T r4 TT i ^ 


and 


i e Ui i e U 


Ci(Ui) = c- -. 2 2 P 

^ D e U, ^3 


i < D 

ijD e 


(2,7.4) 


(2.7.5) 
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Thus the objective function of the current candidate 
problem is given by (2.7.3), ¥e next describe the method to 
determine TESTl, 'JISST2, '£SST3. 

2 

In vie?/ of the fact that y^ = y^ for 5m = 0 or, 1, the 
minimum value of (2,7.3) is same as, 


Min H(UO + I 

i e 


[si-OiCup 


y.- +2 E P . , y . y . 

A y A i d 1 d 

(2.7.6) 


i < 3 
i, j G 


f or y = 0 . or , 1 ¥ i G . 

The follov/ing results are utilised to develop the tests. 


Lemma 2.2 

If for some k G 

_ + Pj, ] >0, 

then y^ = 0 at the optimal solution of the problem (2.7.6), 
Proof ! The objective function can be re?/ritten as 


H(Oi) + yj. 


k' 




+ 2 


2 

G 


f\(k } 


'31c b 


+ 


j G { k } 




b + 2 


E 

i < 


P - • y . y - 
2-0 '^ X '^0 


i,j 6 uA^k} 


(2.7.7) 


Thus , 


^k ” ^ ® 
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=> + 2 


uA{ic> ^ 


^ y^ = 0, or, 1, 

3 e U^{ k } 


Thus the value of (2*7.6) with = 0 is smaller than its value 
with = 1, irrespective of the values of y^ , j 7 ^ k. 


lemjaa_A,*^ 

If for some k g TJ^j 


< 0| 


-c, (U. ) + P, + 2 2 P 

^ . ^ jGTJ^fk) 3k 

then y^^ = 1 at the optimal solution of the problem (2.7.6). 

Proof : Prom (2.7.7), it is clear that 


- Cj^(U^) + P^^ + 2 


D e %\<k 




> "C]^(U^) + Pj^. 


for all y - = 0 or, 1, 
J 


+ 2 2 P y < 0 

0 e u^{k} ^ 

j e u^\{k} . 


==> = 1 at the optimal solution. * 

Thus for the candidate problem under consideration, by 
use of lemma 2.2, 2,3, some more variables coai be fixed at the 
value 0 or, 1, let U^, Uj_, U|. be the set of subscripts of 0-, 
1~, and free variables respectively, after successive use of 
lemma 2 ,2 and 2 ,3 • 
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TBSTl : 



If = 4> , then the candidate problem is fa.tiiomed and 
the optiraal value of the present candidate problem is 


H(Up. 

If 7 ^ i}) , the augmented candidate problem is 


Iain H(U’) - Z Ci(u')y. + I P . yf 
i e ui i G TJi 

i J- 


+ 2 


Z 

i < 3 

is3 e 



y. 


(2.7.8) 


subject to yj_ = 0 or, 1 , ¥ i G U^. 

Since the problem (2.7.3) and (2.7.8) ane of tne same form, 
we shall consider (2.7.3) as the current candidate problem even 
for the case of (2.7.8). 


At- this stage, trivially, 


-c^(U^) + Pj_ < 0 , ¥ i e 

and -c.(U. ) + P. + 2 Z P.. > 0, ¥ i e U^, 

^ 3 e Uf\U> 

becanse, otherwise we can app 13 ^ lemma 2.2, 2.3 .-nccessively, 
till Yie reach such a sta^e or till we fcutaom the candidate 
problem. 


lEiS'1'3 : Select the subscript m G such that 




Min c.(U^). 

3 6%^ 


We shall, use y^ as the branching variable, and two subproblems 
are generated. Ihe candidate problem 1 has the subscripts sets 
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as Uq, {m }, U^{m}, and the candidate problem 2 has 

subscripts sets as U^iJ {m}, U^, U^{ru}. 

For the pruning purpose, the lovjer bound can be determined 
as given belov;, The problem (2,7,3) can be vv-ritten in matrix 
f orm as , 

D(U^) - c(U^)®y + y'^Q(Ui)y (2.7,9) 

vmere , 

c(U^) is the rector (Cj_(U2_)) of appropriate sise, 

y is the vector of (yj_) of appropriate size, and the 
matrix QCU^) is obtained by eliminating ror, and i^^ column 
of Q, for all i 0 U^, 

Clearly, Q(U^), being a principal subuatrix of Q, is 
positive semidef inite . 

The unconstrained minimum value of (2.7.9) can be obtained 
by the method of differentiation, and is given by 


1 


o(Uj^)®Q(up~c(TIi) 


H(%) 


(2.7.10) 


Thus l4(U^) gives the louver bound for the c^timal solution of 
candidate problem (2,7,3). 

The proposed algorithm works as follows? 

Initially, - 
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If the current candidate problem is fathomed then, set 
w = Min (H(U^), w), where li(U^) gives the value of the present 
incumbent and a new candidate problem is picked up from the list. 
If the list is empty the algorithm stops. Ihe subproblem v/hich 
is listed last in the list is chosen as the current candidate 
problem. Other?t7ise, by application of 1TES'12 and 'fSiSlo, the 
current candidate problem is partitioned into tv^o subproblems. 

Ihe subproblem 1 is taken as the candideite problem and subproblem 
2 is listed as unfathomed problems. If LB is the lov/er bound 
of the objective function for the candidate problem, then the 
branch leading from the problem is said to be primed i;hen 
w < 113 . If it is not pruned, then algorithm is repeated v/ith 
the current candidate problem, till the list becomes empty. At 
this stage the value of w gives the optimal value of the original 
problem. 

A. 

Substituting Zj(h) = w, 1 is calculated. 

R emark : Since the proposed branch and bound .algorithm is 

aimed to be used in Algorithm 1, the initial value of v can te 
updated for different values of h from the information available 
about the Zj(h-l). Such improved values of w help to reduces the 
total number of nodes to be scaimed in Algor it nm 2. lue method 
of modification of initial values of w is given below. 

(i) If h = 1, w = (2h+l) 

(ii) If h > 1, w = (2h+l)a'*^-a*'^Qa*-6"’Q&*-Zj(h-l) 
where a* is def ined in (2 .4.2) . 
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Ihottgh in general the branch sna bound 
w 0 ao its initial value, but while solving 
touna problem such modification is very much 


al,.3orithrii assumes 
for Bonferroni 

useful . 


il lustra tion 


Following numerical example 
Al{r,oritlim 2 for n = 6 and 


illustrate 


s Alg or i t iim 1 an d 


0 .62827 

0 .49 564 

0 .20446 

0.49564 

0 .51940 

0.22442 

0 .20446 

0.22442 

0 .29661 

0.31045 

0 .23636 

0.03655 

0.32903 

0.23507 

0.11250 

L 0.09525 

0.0 7645 

0.08274 




Al.SP rithm 1 


0.31045 

0 .32903 

0.09 525 ■“ 

0.23636 

0.2350 7 

0.07645 

0.08655 

0.11250 

0.08274 

0.51358 

0.3500S 

0.08274 

0.3 5008 

0 .57283 

0.00000 

0.08274 

0 .00000 

0.12142 


ii = 1, u = n 


(0.51135, 0.21314 , 0,47680 , 0.44461, 0.75346 , 0.33674) 
I^(h) = 1^(1) = I ^ 9 ^ 

^l(li) 2 !j(l) _ -1.80425, obtained by using algorithm 2 '. 

I' ='0.90213 


k 1. Hence algorithm terminates here, 
liie test Bonferroni Bound is 0.90213 at 

= 0, ag = 1, ^ ^ 1, h = l 

A lgorithm 2 

Initial lo^/ver bound is -I.91I33 
Initial upper bound = w = .- 1 . 69896. 
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= { 1,5 } , Uq = <{. , = I \{1,5} 

v^hicb has the value LB = -1.84714, and the nev.' entry to list is 
the candidate problem 

= {1} , = { 5 }, Uj = n\{1,5} 

which has the IB values as -1.77521. 

Using OZFjSII, IE,312 successively, the current cendide/te problem 
is fathomed. Tne complete solution is 

Yl = y2 " ^3 "'O’ ^4 "" 0’ ^5 " ^6 " ^ ! 

and the value of the objective function s.t this point is -1.79647, 

So the value of w is augmented to 
w = -1.79647. 

hith this new value of w the list is pruned and latest candidate 
problem is deleted from the list, 

Thus the algorithm is repeated till the list is empty, The 
whole algorithm is illustrated in fig, 2,1, 

2. 8 Oomputational vesutts 

In this section computational experience for Algorithm 1 
and Algorithm 2 is reported. Ihe algorithms are coded in 
]?0RIRAH-4 and run on lEG-lO (KL4o Processor) system. Algorithms 
are tested on randomly generated problems. 

Pirst we report the computational results for the branch 
and bound algorithm. In table 2.1, the results refer to average 
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LB = _i .77521 
Bruaed 

i I'B= / 

-1.84714 \ / 


fatbomed 

conplete solution 

=1,y =1 


w = -1 .79641 


- 1.91133 
-1 .69896 






= -1.88011 


= 1 

/ 

/ 


^5 = 0 



\ LB=_i,67416 

\ Pruned 


)lb=_ 1.876 13 / 


^2 = 1 



LB=~i .83712 



O 


LB = _i .76027 
Pruned 


fathomed 

conplete solurioa 
^1=0,72='’ »y3=0,y^=0,y^=1 ,yg=i 
W = -1 .80425 


Optimal value = -1.80425. 

Time takeaby DEG System lo = o.08 aec. (C.B.u.) 

2.1 ; BRAUCH Affl) -Rnnivm 


49 


time and average number of nodes scanned, for randomlj generated 
problems of different sizes as indicated. 


Value of 
n 

No. of problems 
s olved 

No. of nodes 
encountered 

lime talren 
( 0 .P . 0 ) in se 

6 

50 

7 

0 .06 

7 

50 

8 

CO 

o 

• 

8 

50 

8 

0.09 

9 

50 

12 

0 .12 

10 

30 

15 

0 il6 

20 

30 

40 

0 148 

30 

20 

50 

1.32 

40 

20 

77 

5 . 54 

50 

20 

400 

30 

loo 

10 

1263 

135 


Table__2 _._l 

Next vre discuss the con^utational experience for 

Algorithm 1, It is observed that for all the problems generated 

randomly, the approximate value of L obtained by Algorithm 1, 

is better than Max L- and thus for such problems Algorithm 

1 < i < 8 ^ 

1 gives the best lover bound. 

In table 2,2 we give some typical problems of different 
dimensions. Ihe average number of iterations and the average 
time taAen by the algorithm is also given against the respective 
values of n. 
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Yalue 

No. of 

Average 

Sirae taicen 

POE A SIPICAl PEOBiEii 

of 

problems 

number of 

(0.P.U-) 



n 

solved 

steps 

in se cs , 

Approx# 1 

L* ! 

. 1 

-1-1 » 

1 < i < 8 ^ 

6 

20 

3 

0.16 

0.94145 

0.92753 

7 

20 

3 

0 .20 

0 .8854:6 

0.84317 

8 

20 

3 

0.31 

0 .850 6 5 

0 .82429 

9 

20 

3 

0 .40 

0 .89357 

0.89342 

lO 

10 

5 

1,01 

0.95941 

0 .9 5509 

20 

5 

7 

3.57 

0 .95544 

0 .90859 

30 

2 

15 

90 

0 .9 3904 

0.89953 

40 

2 

20 

130 

0 .90009 

0 .90009 

1 s 
i 

2 

20 

180 

0 .88608 

0 .88608 

Sable 2,2 


In table 2«3, we give for different problems for n = 6, 

the values L and Max L . . I'ne correso ending Q matrix is 

1 < i < 8 ^ 

given in Appendix !♦ 


Problem 

number 

Approx. li* 
by Algorithm 1 

Max ii . 

1 < i < 8 ^ 

1 

0 .9 76 

0 .962 

2 

0 .555 

0.54833 

3 

0.72 

0 .695 

4 

0.88466 

0.85147 


Sable 2,3 
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dravmj 


From the above results two taportart oonolusions osui be 


(i) Previously there was no systenatio method to compute a 

r bound lor high values of n. The present algorithm provides 

an efficient method to compute the bound for reasonable values 
of n. 

( ) looking at the values of the improved boim.ls and at the 

n»ber of steps taken by Algorithm 1 , it is believed that 

justified by approximating L*. 

2.9 Computation of lower hound L 

B 

It ie seen from corollary that .ha:: ,130 gives 

a ound tetter than the lov,er bounds obtained m [; 31 

this section it is shown that oo^utation of such test’ ilovn: 

Wr bound leads to n :^kimisation problems of ocncave „io 
Pseudo-Boolean function of 2n variables. 

Su.l}sij in^ ( sl f' nY' Y" t 

S ^a,n; lor Z, Lg can be determined as| 

,* n 

- Max 2 p X - 2 n 

Z e y* i=i i i 1 < i< j < n ''=13 


-1-.X 

(a,h) e yI 


a’^ - a‘^Qa 
h(ii+'l) 


( 2 . 9 . 1 ) 


An in (1.5.1) . (i.s.yj, ,,3 ^ restricted to 

1 i I- < n, from oorollary P.l and 2.2, it is clear that we need 
only to find the quantity 

: *- ■ V (. I, t . 


Max 

1 < h < n ® ^ 







which imprcves the hounds 1^, 2 < i < 7. 

Putting 

ai = 2y^ + - 1, 1 < i < n (2.9.2) 


v\?here, yj_ = 0, or, 1, = 0, or, 1, the problem of de terming 

Lg(h) reduces to 

Maximize H¥+TT I (2h+l)(2y+z-^ )'-^-(2y+z--e )'^Q(2y+z-e ) 
y e z e Bg ^ 

This 0-1 maximization problem is equivalent to 


P3 : 

where 


muiirnize I E 

a e Bg 

, e is a vector havit 


Qe + ( 211 + 1 ) 0 "^ —> 
coirp onents 1, 


T 

f a + a 





f 


T 


(2(2h+l)P + 20^0)^^, ((2h+l)P + e^Q)'^‘ 



P •r)2n 
8 Bg , 


H = 

It is easy to see that, 

Q is positive semidefinite ==> M is positiTO semidef inite . 
Hence the objective function of P3 is convex, iloreover, M ■ has 
nonnsgative of f -. 4 iagonal elements. * 


4:Q 2Q 

,,,2Q Q J. 



Hen<» the problem P3 can also be solved by Algori% 2 



2,10 Some properties of v* 

JR 

It is observed from r? p -v, + 4 ., 

(2.2.2) oiiat the integer restriction 

rise to the class V*. UnliVP +h * 

“ase of 7*, any x e V? need 

not aato^atioally tdat xe P(.). for Instance. 

n = 2, h=i p, 

» ±, - 1 + ag = 1 + Soj 

where e^, Sg are in open interval (n i ^ +- 

^ interval (0,1), the point (a,h) e Y* 

does not satisfy all the constraints (2.2.a). On che contrary, 

If Bp eg lie in the close-open inteiYal the associated 

point (a,h) s is feasible, 

Thus clearly, 

Y*ih) C Y*(h) n P(H) 


==> l*(h) < 


X e n < g ^ j. 


h is fiiied. 


j"(2h+l)a^ - 
(a,h) e VgCh) n P(ii) L'’™™^F(h+~ 


=> Ii*(h) < Max I 


1 . 

a ya 


( 2 , 10 . 1 ) 


b is fixed 


< 


Ihus L* < Max l^(h) 
“ h > 1 ^ 


lemma 2.4 gives the value of Max L^^h) 

n ^ 1. . 


( 2 . 10 , 2 ) 


(2.10.3) 
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Lemma 2.4 


Max Lr(h) 
h > 1 ^ 


9 T - 
I P Q P 


Pr oof : By the method of differentiation, we have 


(2h+l)2 _ 

L (h) = P^QP, 

4h(h+l) 

and it is attained at (a = Q~P,h). 

Since Q is positive semidef inite , so is Q“. Hence p'^^q'P > 0 
Aga-in being a decreasing function of h, for h > 0, 

( 2'ri+iO ^ _!E_~ . . - _o j. • ■ 

/i"7T'’tT P Q P IS also a decreasing function of h, 

4n(n+l) 


Hence ^Max^ =1^(1) = § P^Q~P 


phis completes the proof 


Thus from (2,10.3). 


< I 

The following observations are apparent, 
3 


(2.10,4) 


(i) If ^ Q P is an Integer vector then 

l* = |pVp. 


This is so because the point Q P,l) e where L^(h) take 


s 


9 ^T 


its maximum value ^P QP. Hence, from (2.10.4) vje have 


1* = § p®qT> . , 


(ii) If for a=^QP, h=l, (a,h) 6 is a fea,sible point of 
P(H), then the best known lower bound is P Q P , 
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This is directly, implied by (2.10.1), (2.10.3). 

Ihus lemma 2,4 gives a lower bound of 1*. ‘Ineoreiii 2.7 
is useful in finding an upper bound of L*. Besides this, it 
gives an insight to the geometrical structure of Y*. 

ll 

Theorem 2.7 

lor any fixed b e the set of points 


, 2h+l 


h > 0 


n(n+l) 

2 

is a line in R . 

Pr oof : We prove the theorem by shov/ing that the point 

2ii ' -f 1 -jf 

( — b,h ) 6 lies on the line joining the origin and 
(^iLni b,h, h’-l) G Y^. Thus the set of points 

i,h) i h > 01 Cy|, 


gives rise to ray emanating from the origin. 


let x' = tjh'), 

x''= (Sriiu 1, h'.), h''= h' - 1 . 

We claim that, 

z’ = X x" , 


Y^here X 


(2h'+l) ^(h*-l) 
(2h'-l)2(h’+l) ’ 


To see this, 
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x" = 

(gii'-i)(gh'-i) ^ 

. 1,2“ 
4 r_ 

1 

i 1 

L 2 

i 

h«'(r*“-r) 

(2h’~l)2 

T' 


1 


xxn = 

(2h»+l)2 

2h’(ii'+l) 

b? 

= xl 

1 


Similarly, X . = x!.. 

X J X J 




It is not difficult to check that 


0 < ^ < !• 


Thus for any arhitrany h’, 


X* lines on the line joining x” and the origin. 

This completes the proof. 


Corollary 2 .4 

For h’ > h" , C![;v^(h‘)I3C G[[V*(h”)3. 

P roof : Since o G ¥^(h), ¥ h, any point (a,h*) G is the 

convex combination of 

(fl™ a,li") e 7|(h'') and tHe carigin. 

Corollary 2 >5 

For any h > 1, if Q”p,h) e F(H) then 

Q-Pjb’) e F(H), ¥ h* > h. 
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J Since F(ir) is a con-rox set and the origin is also in 
P(N), by using the result of theorem 2,7, the result follows. 

Prom the aboi/e results it can be said that given an 
arbitrary vector a e R^, there alvrays exists a h(a) such that 
for any h > h(a), the point a,h) e P(N) and for any 

h < h(a), the point a,h) gf P(N). 

let h denote the value Q'p(h). Then clearly l^(h) 
provides a lower bound. 

We fail to use the above results as a tool for easy 
solution procedure, Nevertheless this premises a lower bound 
Ig^(h), It is believed that a feasible direction search 
technique can possibly be developed making use of the above 

' ■ 'Ilf 

results. 

2,22 Coham, g&nevabion techniqa to solve the dual 

In this section it is shown that the use of column 
generation technique to solve dual of IPCN) leads to maximizing 
a sequence of Pseudo-Boolean quadratic functions* But unlike 
the cases of P2 and P3, such problems can not be solved by 
Algorithm 2 discussed in 2,6, because the method assumes the 
function to be concave with positive off-diagonal elements of 
the associated quadratic form, which may not be true in the 
case of the subproblem generated by column generation* 

Writing the primal linear program IP(N) in the matrix 
form, we have 
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Max P^I 

LP(II) : 


( 2 . 11 . 1 ) 


subject IK < e 

where e is a ■vector ?/ith all components 1, P is the coefficient 
matrix whose rows are 1 ( 0 ), ¥ ^ ^ i'* • c\sGoci;'ted dual 


linear program IP (IT) is 


. . 1 

Minimize e I 

IP(N) : m 

such that D Y = P 


Y > 0 


( 2 . 11 . 2 ) 

(2.11.3) 


Since corresponding to every subset JC K, there is a component 
of the. dual variable -vector Y, we denote the components of Y as 

y(J), ^ J SK. 

Ihe column generation technique stsu^ts ■.■.‘ith a known basic 
feasible solution Y^ of IP(H). let 3 be the associated basis. 

I hen cle arly 

Yj = P > 0 (2.11.4) 

Prom the theory of linear program Yg is optimal iff 

e® 3 “^ d(J) -1< 0 (2.11.5) 

for all J C IT. 

If Yg is not optimal i.e., if Yg does not saxisfy (2.11.5), then 
a column d(J) of D is chosen such that 

e^ B-^ d(J) - 1 > 0 


( 2 . 11 . 6 ) 
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Uhen the value of the objective function of 110(11) can be 
increased by introducing Y(J), satisfying (2.11.6), into the 
basis. In view of the fact tha,t the coefficient natrir contains 
2^-1 columns, the problem of checking the criteria (2.11.5) and 
(2.11.6) les.ds to a formidable task. Column generation techniq.ue 
handles such problem by formulating as a suboptimization problem 
as follows? 

I -1 

Maximize e B Z - 1 

SP(II) : 

subject to Z e {d(J) j 4 ) ^ J C il}. 

The suboptimization problem SP(H) is eq.uivalent to a 
linear ‘■progr am where the aim is to maximize the linear function 
over the convex set 0 'which is the convex hull of the set of 
points ^d(J) I ^ ^ 3 

Theorem 2.8 sho?;s tnat such a convex set is a convex 
polytope of diameter 1. 

Theorem 2.8 

Every d(J), for J ^ is an extreme point of 0 and every 
pair d(J‘), d(J”), J’ J” are adjacent extreme points in C. 

Proof : The proof given here is based on the similar type of 

pr oof given in ^ H • 

let for any J,. |J j = k, 1 < k < n, toe vector 
n(n-fl) 

f(J) e R ^ be defined as 
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such that = 1 if i e J, 

= 0 otherwise , 

= -1 if i,3 e J 

= LI > 0, otherwise with i: > 0 OiT'bitrc'Xj.u 
Clearly the inequality 

f(J)^X < holds for e-very Z e { d(3); 4) S C_ N} 

and equality holds only if 1 = d(J). Inus, 


1 

f(J) X = 


k(lc+l) 

”'2™' 


is a supporting hyperplane to C v^hich intersects it in d(J) and 
hence d(J) is an extreme point of C. 

Next we shall shov/ the adjacency among any pair of extreme 

points. 

Nor a given pair of subsets J', J ”, J’ J ” let us 
de f ine t he f oil ow ing , 

s ' = J ’ n J ” , s ” = J ' ■ y ” , s = J ” \j ' 


with is* j = s^, |S”| = Sg, is*” 1 = Sg, jJ* j = We determine 

n(n+l ) 

the vector f(J*,J’‘) 6 R ^ as follov/s, 


f ( J * ,J '*) • ’^ln’^23* * * *’^n-ln^ 


where , 

= 1 if i e J*, 
^2 

= -^ if i e S ”* , 
®3 

= 0 other v^ise , 




61 


and f.. = -1 if i,3 G J' 

X J 

= - if j e s"s i G S’ 

SgCSg-l) 

- - sTCig-TJ e S’” 

o 

= -rii othervvise, v/itli M > Ic^. 

Let us consider the inequality 

r|i 1- "1 ( 1- 'I 1 ) 

f(J’ ,1")-^ X < — 

2 

Y>?hich holds for all X G ■td(S) : 4i / S ^ II } ajid e(|uality holds 
only if X = d(J’) or, X = d(J”). ilence the hyperplane 
X^(k^+1) 

f(J‘,J”) X = "---—““ intersects 0 at two points namely d(J' ), 
d(J”)* Ihus d(J’), d(J”) are adjacent in. C. 

this completes the proof. 

Ihus it is not fruitful to solve suen a suboptimization 
problem by usual linear programming technique. I'liis is so 
because the number of extrene points that are adjacent to any 
extreme point is exponent tally large and any method trying to 
select the best among the adjacent points vi?ill require o(2^) 
number of computations. 

next ?;e shoY\7 that SP(II) can be considered as maxim tza.t ion 
of a quadratic Pseudo-Boolean function. 


let 
I- —1 

e B = (q^jig? • • • 


Si-ln) e B 


n(nj-l ) 
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I>efine u e f. 


X e {d(J) I JO at} ^ ^ 

-^±i clearly 

- ^ir SP(R) is equivalent to 


Max _ ^ 


Spl(H) : 


1 / - t . li a . u . 

^ 5, 1 < 3 < n ^ J 


“l = 0 or, 1, 1 < 3, ^ 


Jo special purpose technique to solve SPl(H) is 

discussed here, in order to solve the tv c i 

technique, it is neoulr- i + *' ^ ol»dh fjeneration 

, ®1C0 se^teral times. Por 

large values of n s„ct, • i or 

. ^uch an ^proach turns out to inefficient. 




Chapter III 


b:ithstb__ p pjiras_, _pp 

Inti?oduat-ion 

Study of the geometrical structure of tue constraint 
polyhedron of a linear program is a -yery valuable information 
for de ye loping its solution techniques. In generaly the 
geometrical study of the polyhedron concerns nlth the structure 
of its extreme points and of their adjacency. Such an invest i— 
gelion is more important, if the linear' program has a, large 
number of extreme points. Keeping this in mind, study of the 
structure of P(K) is initiated in this cnapter. Ihough we fail 
to supply the results of this chapter in coiri^utin^ the lower 
hound, we stress on the theoretical potentiality of these results 

Ihis chapter mainly deals v/ith the geometry of the 
convex polyhedron P(Ei). In particular, we explore the structure 
of its extreme points, first, one type of extreme points are 
identified. It is shown that such extreme points are memhers 
of V*. Ihe necessary and sufficient condition that any point of 
7* is an extreme point of f(N), is also derived, Anocher type 
of extreme points of f(l) is characterized and it is shown that 
such points can he obtained from certain nonextreme point 
members of 7*. «'e believe that these t\/o classes of extreiae 

points characterise all the extreme points of f(U), but mq are 
unable to prove it. Attea^ts mads in this direction are reported 



and in the process the class of Kounias and I'ariii, given in 
(1,5,1), is generalised to a broader class. 


The adjacency structure of soiiie of the hrio-vm extreme points, 
is also discussed in this chapter. It is shov.n that the points 


8 

no't in the above two classes 
i=l 

are not extreme points of 3?(li)* 


of extreme points, 


We next introduce some preliminary concepts for immediate ; 
use in the subsequent sections. 

2.2 PritimCnary oonoepts i 

It is seen from (2,1.2) that there is a one— one correspon- 
dence between the family of non-empty subsets of II = {l,2,...,n} ' 

and the set of constraints of liP(N). Thus without any ambiguity, I 

v?e shall identify aiiy constraint by the associated subset of N, 

Let.S be a subset of N having s elements, We assume without 
loss of generality, that S = {l,2,,,.,s> . Let T, J, I denote 
any subset of S,1I, s\(i' respectively. The subscript attached to 
T, J, I ?/ill denote the cardinality of the subset. I'or instanc®, 

T^ denotes a subset of T of S for which jlj = k. Such a 
subscript is used ?/hen the cardinality of the set is needed for 
discussion and dropped when it is irrelevent, G-enerally Greek 
letters will be used for fan Hy of sets. 


Let 6jj.(S) denote the family of all subsets of S having k 
elements, for a given I Cll\S, let 
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Examp le 5«1 : 

Let 1 = {1,2, 3, 4,5,6 } ,S ={1,2,3}. 

Let I = {4} , then I'^(S) is given by i 

{4}^(S) ={ {1,2,4} , {1,3,4 } , {2,3,4} . 

Similarly when 1 = {4,5}, 

{4,5}^(S) ={{1,2, 4, 5}, {1,3, 4, 5} , {2,3,4,5} . 

^ Jlow the set of constraints (2.1.2) is arranged in a 

rectangular array R(S) consisting of (s+1) rov7S and of 2^*“® 
columns. Each column corresponds to a subset I C L\S and the 
rows of R(S) will be counted as 0,1,2, ...,s. i'ne entry in 
row and colunm of R(S) contains all the constraints associated 
?7ith the block I^(S). 'thus each "entry of R(S) contains numter 
of constraints of (2,1.2). It is to be noted that 4>°(S) is 
empty. 


Example 3,2 : 


Let H, S be defined as in Example 3.1. fnen the 
rectangular array R(S) is given by 


k 

0 


<}> 


(1)(2)(3) 


(12) (13) 
(23) 


(123) 



:5___ _46_ 
(45) 1(46) 


(145) 

(245) 

(345) 


1245 

1345 

2345 


i 



(146) 

(246) 

( 346 ) 

1246) 

1346) 

2346) 


. 56 ^. 

(56) 


(156) 

(256) 

(356) 

(1256) 

(1356) 

(2356) 


456 


(456) 


( 1456 ) 

(2456) 

(3456) 


12456 
; 23456 
' 13456 


(12351 (1236) kl2345)j(12346)| (12356)! (123456 


1 •. 


^ 5'iil 
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As in section 1.3, let 

n(n+l) 

■mm : . nwXm a iiit Ji mAm 

d(J) - ‘••’'^’‘^12’ ‘••»%n^‘^23’ •••’^-In^ ® ^ » 

be a vector defined for each nonempty subset J of 1. Given a 
J i^be first n components of the vector d(u), Vi’ill be denoted 

by the vector 


n 


d(J) = (d^jdg, e E 

definiti on 5*1 ; Ranks of a family of subsets 

Eor a collection a(ll) of subsets of N, v;e define rank of 
ord er two (ran k of ord er o ne) as the maximum number of independent 
vectors in the set 

{d(J)jJe a(N)> ({d(J)jJe a(H)>) 

and will be denoted by r 2 (a(l)) (r^(a(K)) ) , 


Clearly rg(a(N)) < , and r 2 _(a(li)) = n. It is to be 

noted that whenever the collection a(lT\^S) is referred the 
constraints under consideration are that of IiP(E\s) i.e., the 
linear pr ogramning problem with respect to the events A^^, 
i G e\s . Similarly in this case, the vectors d(I), d(I) are of 
dimensions , (n-s) respectively. 

lo see the relation between R(S), a(l\s) and the set of 
extreme points of E(N) let us consign the foliating examples. 

Example 5.5 ; #■’ 

N = {1,2,3>, S = {1,2> • Ihen E(S) has 3 rov^s and 


2 columns 
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0 

1 

2 




an (2) 
( 12 ) 


) I 7 




^2 

(3) 

_(i2)/(23y 

(123) 


For n = 3. thprp a-ro -P 

^ , , . _ extreme points of ^(11) and all such 

erne points are obtained from the folio, i nr- v-ir nf - ■ 

u-i-Lu.,ina Juir of pomrs 

y permuting their indices ^ 31 2 ]. 


1 


Xr 




1 

1 


^12 

1 


-^3 


^23 


^ 1 1 1 

(3.1.2) 

° ^ 0 0 ( 3 . 1 . 3 ) 

I'he following ane set of constraints that are satisfied as 
equalities at the point (3,1,2) 

d(S), ths), {31°(S), {3h(s) (3.1.4) 

Similarly the constraints which are equalities at (3.1.3) (rj 

♦ts), /(S), (31(3), {3)®(S) (3.1.5) 

Tor other extreme points such set of constraints can he written 
by permuting the indices of » in (3.I.4) and (3.1.5). 

let E = {1,2,3,4} . For n = 4 all the extreme points 

can he found from the f dicing 5 points by permuting their 
indices . 
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^1 

Xg 

^3 

^4 

^12 

Xfg 

^14 

^23 

""24 

^34 


1 

,1 

0 

0 

1 

0 

0 

0 

,0 

0 

(3.1.6) 

1 

' 1 

1 

0 

1 

1 

0 

1 

0 

0 , 

(3.1.7) 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

(3a.8) 

1 

1 

1 

-2 

1 

1 

.1 

I 

-1 

-1 

(3.1.9) 

2 

2 

2 

2 

1 

1 

1 

1 

1 

1 


3 

3 

3 

3 

3 

3 


3 

‘3 

3 

( 3 . 1 . lO ) 

Por S 

= {1, 

2> 

the 

associated R(S) 

' is 

wr it ten 

r-" rn 
Ci,0 

follows 

1 



(3) 


(3,4) 1 

Cl)(2) 

90 

^ (23) 

(14J 

(24) 

( 134) I 

(234) 

'“'cTsl 

(123) 1 

J 

(124) 

(1234) 


In terms of the entries of R(S) the constraints which are 
equalities at the point (3.1.6) - (3.1.9), are respectively 

<(,^(S),<j>^(S),{3 >^(S),{3>^(S),{4>^(S),{4:}2(S),{3,4>^(S),{3,4}2(S) 
<()^(S),({.^(S),{3>°(S),{3>^(S),{4>^(S),{4>^(S), {3 ,4>°(S ) , { 3,4 }^(S) 
(j)^(S ) , <}.^( S ) , {3 >°(S ) , ( 3>^(S ) ,{ 3 >°(S ) , {4 f (S ) ,{ 3 ,4 >°(S ) , {4>°CS ) 
<i)^(S),/(S),{3f(S),{3>^(S),{4>^(S),{3,4>^(S),{3,4>2(S). 

In a similac manner if S is taken to he H then the point 
(3.1.10) can he obtained hy making the constraints 6 (W), (N) 

as equalities. 

Prom the above examples, it is observed that given an 
extreme point of P(H), n < 4, a suitable S can be chosen such that 
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one associa,ted E.(S) has the following properties* 

( i) Ihe extrerae point satisfies either all the constraints 
of a particular block t"(S) as eciualities or, none of the 
constraints of the block is satisfied as eauality at the given 
p oint , 

( ii) Eor a given 1 q. i.e,, for a given colurm of R(S) 

the extreme point satisfies at most tv/o, blocks as equalities. 

In other vv^ords, there are at most two integers , ko such that 

Ir V ^ d 

Up 

the blocks I (S) and I (S) are equalities at the given extreme 
point, where O < k^^ <_ s, i = 1,2,, Further, 

(a) if the point satisfies exactly two blocks of column I, 
then the blocks are consecutive i.e., 

jki*-kg j =1. 

and (b)if the point satisfies only one block of I then tte block 
is either the first block or, the last block of the 
column of R(S) i.fe., either = o, or, k^^ = s. 

( iii) ‘Ihere exists a collection a(N\^S) of subsets of such 
that 

rg(a(n\s)) = , 

and for each I e a(N\^S) the associated column in R(S) has 
exactly two consecutive blocks as equalities for the point 
under consideration. 

In the next section W show that for any arbitrary n, 
we can obtain a class of extreme points of F(W), by suitably 
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ohoosing S and following the rules (i), (xi) aiid (ill) which, 
in fact, have been observed for the case n < 4. 

3,3 ExtCl) : Extreme 'points of first hind 

1‘rom definition 1.1, it is clear that ai'd' entieme point 
of I'(IT) is a solution to a system of linetiT e neat ions obtained 
by considering some independent const-' c.i-rn. of (2.1.2) 

as equalities, i'he method adopted here to generate extreme 
points is to select a set of constraints from (2.1.2) such that 
this set contains a set of linearly independent cons- 

traints as its subset and that it follows the rules (i), (xi) 
and (iii) for some given S. This set of constraints when 
considered as equalities gives rise to a feasible solution 
(obviously, unique) and hence, this solution is an extreme 
p oint of l'(h) . 

This section contains 5 subsections. In the first three, 
we gradually develop thi| results which are useful to prove the 

if’ 

main result of this chapter. The main theorem is given in 
3.34. In 3.35, we generalise this result and e,ive an 
illustration. • 


3.31 In view of the rules (i) and (ii), here we consider the 


case Vi/hen two consecutive blocks of column (!>(S) of Pl(S) are 

• equalities and derive a system which is equivalent to it. 

' : 

let, for some integer h, 1 < h <_ s-2, (!j^^(S), (>^"^^(3) 


be equalities. The system 't)f equations so obtained is. 


X 

i e T, 


X < D 




(3.3.1) 


I','., / . 
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for all e 6i^(S), and 


r x- - Z 

e 1 


Xi- =1 


(3.3.2) 


f or all 1 


h+l 


e 6 


h+l 


h+l 


(S). 


If is easily seen thaf system (3.3,1) occurs for 1 < h < s 
and system (3,3.2) occurs for 0 £ h < s~l. The comhin*ed system 
(3.3,1) - (3.3,2) occurs for 0 < h < s. 

Theorem 3 ,1 : 

The system of linear equations (3,3,1) — (3,3,2) has the 
following point as the unique solution, for all values of h, 

1 < h < s-2| 


h+l ’ 


X. . = 


^3 h(h+l) 


> i <Dy i, D G s 


(3.3.3) 


Qo t i By direct substitution of values of x-, x. . from 

i i 3 

(3,3.3) in (3,3.1) - (3,3,2), it is clearly seen that (3,3,3) 
is a solution of (3.3,1) - (3,3,2), It remains to show that 
this solution is unique^ That is to show that among the total 
number of equations of the system there are 
linearly independent equations. In theorem 4=, 7 of Chapter IV, 
existence of such a^linearly Independent subsystem of (3,3,1) - 
(3,3,2) is shown. Bor direct reference, an equivalent statement 
of theorem 4,7 is given here, by substituting II by S and J by T; 
'' There exists a collection ‘a(S) of subsets of S such that 


ot(S)c=6i,(S),__, 


aiiii 
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rg(a(S)) = , for 1 < li < s-2’' . 

jDhis completes tlie proof. 

G orollary 3.1 ; 

For h = 0 and h = s-1, although (3,3.3) is a solution 
of (3,3,1) - (3,3,2)j the solution is not unici_ue , 

• Since |6jS) U6i(S) j = i6i(S)j = = s. Hence 

for h = 0 the system has s linear equations in variables, 

so its solution cannot be unique. Similarly, 

|6s-l(S)U^s(S)| = \_i + \ = s + 1, 

SO the system cannot have unique solution. 


Prom theorem 5.1, it is clear that the system (3,3.1) - 
(3*3,2) is equivalent to the system (3,3,3). 


5 ,52 Here we derive an equivalent system for the system of 
equations obtained when two consecutive blocks of any arbitrary 
column of R(S), considered as equalities, together with ^^(S), 

n(n+l) 

2 

let us consider those points X G R which besides 

being solution of (3,3.3) satisfy exactly two ccaisecutive blocks 
of a particular column I of R(S), as equalities, let the blocks 
which are equalities at the point under consideration be I®(S) 
and I®''"^(S) for some 0 < 6 < s-1. Substituting the values of 
X., X.., i < 3, i ,3 e S from (3,3.3) we have, 

i Ij 


2e . 


i e I 


^i 


I 

i e 




13 


3 G I 


- 2 
i < 3 

3 6 I 


e(e-i) 
^ij : HhTij 


= 1 , (3.3.4) 
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for all Tg e 6 q(S). 
S imilarly, 


^e+ii 


+ 


i e I 


2 X- . 

i e 1, ^ 


" 0+1 


3 e I 


for all e 6«^i(S). 


e+1' 


E 

i < 3 

ij3 e I 




= 1 


(3.3.5) 


Lemma 3.1 : 


The system of equations (3.3.4) aiid (3.3.5) is equivalent 
to the system of equations (3.3.4) and (3.3.6), where (3.3.6) is 

J g J +0 = Iffel) > all 1 e S. (3.3.6) 


l££2£ • ^0^ any Tq^ 3 _ select Tq<^ and let k e 'i'0+^^0 • 

Then (3.3.5) for each Tg_|_^ can he written as. 



E 

i e I 



1 e I I 

3 e I 


E X. . = 1 

i < 3 

i,3 e I 

Subtracting (3.3.7) from (3,3.4) we get, 


(3.3.7) 


36 I 


""kj " H 


(h^) 

traj 


Such an equivalence relation can be obtained for each of 
the constraints in (3.3.5) and hence (3.3.4) ~ (3.3.5) is 
equivalent to the system (3.3.4) tc^ther with 



This coiDpletes the proof, 

5 • Here we derive 

conditions sucn tbat any block E(S) can 

Toe rep)re sBrte h Tov o, sin/yin 4, • » 

S constraint. Once toe blocks are 

reduced to ctogie conetrainte, it le aleo shewn that every 
ooluMi Of E(S) has a donh.aht constraint in the following sense. 
DefmHlon 3^ : Dominant constratot. 


A 

dominant 

(i) 


(ii) 


oonstraint in a set of constraints is said to be a 
constraint of the set, if 

^is constraint is satisfied then the rest of the 
constraints are automatically satisfied and 
a point satisfies any constraint of the set as 


equality, then the dominant constraint is also equality 
at the point and otherwise, the point is not feasible 
to the set of constraints. 

17e now prowe lemma 3.2 which gr^s the condition that any block 
is reduced to a single constraint. 


lemma 3.2 ; 


Let a(l\s) be a collection of subsets of h\^S such that 
r^(a(N\s)) = n-s (3.3.i 

and let tliere be integers h, a^, i e n\s such that 

0 1 ^ ^ ^ I aq < s-l, for all I e 

Then each of the blocks of E(S) is reduced to a single 
constraint if the constramts , ■ 
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^\S), <^^^■^^(3), (3.3.9) 

are made equalities for all I G a(ll\3) where 

0(1) = h - Z a. (3.3.10) 

i e I 

P roof ; By making the constraints I®^^^(S)j I®^'^^'^^(S) 
equalities for I e a(K\^) together with , 5 "’'^^(S), the 

system of equations so obtained is equivalent ’uo, 

(3.3.3) - (3.3.4), for all I G a(l\S) 

\ 

(3.3,6), for all I G a(r\s). 


Prom (3.3.6), we have 


Z , ¥ I G a(K''v^) and ¥ i 6 3. 

i G 1 

Since rq(a(]i\^)) = n-s, the integers a^, i G are uniquely 
determined by (3.3.10). Hence v«/e have. 


x-y,. \ 

=^ij = K(mT) » ^ 3 e and ¥ i G 3. (3.3.11) 

How let us consider any block 1^(3) for any By 

substituting the values x.,x-^ from (3,3.3) and (3.3.11), it is 

'V.; ■ ^ ^ . 

easily seen that the whole set of inequalities in 1^(3) reduces 


to a single inequality given by. 


h+T ^ q g h(h+l* 


- k Z 


i G I 


2a^ 

h'(h+l)' 


Z X. . < 1. 
i < 3 ^ 

ij j 6 I ( % :■ 


(3,3,13) 


Shis completes the proof 
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In theorem 3.2 we show the existence of the dominant 
constraints . 

Theorem 5,2 : 

The set of constraints associated v/ith any column I of 
R(S) has a single domjinant constraint under the conuition 
(3.3.3) and (3.3.11). 

proof j Prom lemma 3.1, the condition (3,3.3) and (3.3.11) 
reduce any block to a single constraint and this constraint 
with reference to block l“(S)| is given by, from (3. 3. 12); 


2 X. 

i e I 


2 

i < 3 
e I 


^13 




m + Htraj” i g j Errar 


Let Eii(m) denote the right hand side of (3.3.13) i.e. 

RH(m) = 1 - + m __ 2 ^ E^rnTJ * 


i e I 


T hen we have , 

RH(m) ~ RH(m+l) 


2 2m V 

E+T - ETE+TT *■ i I j EIIhTJ 


(3.3.13) 


Thus, 


hi'Evry 


( h ~ m - 2 a - ) 

i e I 


(i) for m < h - 2 a., RH(m) > RH(m+l) 

i e I 

(ii) for m > h - 2 a^, RH(m) < RI-I(m+l) 

i e I ^ , 

(iii) for m = h - 2 a-, SS^a) = RH(m+l) . 

i e I ' ■ 


(3.3.14) 
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Thus RH(m) taJses its minimum value at 

m = 0 ( I ) = b - 2 a • , and atm = 0 (I)+l, 

i e I ^ 

For a fixed I, left hand side of (3,3.11) remains constant for 
different values of m and hence we haAre folloninc two cases , 

Qg:§.g_i.A l J If 0 < h - 2 a- < s-lj then clearly the 

i e I ' 

inequality (3.3.13) corresponding to m = e(I) or, e(I) + 1 is 
dominating among I^(S), 1 < m < s. 

Case (ii) t 

(a) If h - 2 a- > s-1, then clearly the last block 

i e I 

of the column I(S) is dominating over all the 
constraints of the column, Tnat is I^(S) is the 
dominant constraint, 

(b) If h ~ 2 a. < 1, using same argument I°(S} is the 

i G I 

dominant column of the set of constraints - I(S)* 

This completes the proof. 

Remark ; In theorem 3.2, it is shown tha.t under the condition 
(3,3.3) and (3,3.11), the dominant constraint of a column appears 
in exactly one block for some column and for others it appears 
in two consecutive blocks. This result helps to consider tte 
rule (ii) mentioned in section 3.2, for generation of extreme 
p o.ints . 

3 .34= In this section we gi-^ the main result of the chapter. 

By using the results obtained in sections 3,31 - 5.33, v/e derive 
conditions to generate adass of extreme points. 



78 


Since all tlie results of sections 3.31 - 3.33 are obtained 
for a specified S, there is no loss of generality by dropping 
the suffix S from we henceforth, use instead 

of 

Theorem 3.5 : » 

Let the integers h > 1 and a^, i e be chosen such that 
they satisfy the follovving conditions^ 

( i) 1 < h < sr-2 

(ii) there exists a collection of subsets, a(N\^S) such that 

rgCaClN^)) = )a(N\s)j = 

( iii) 'o < e(I) < s-1, 

¥ I e a(N\s) where e(I) = h - 2 a.. 

' i e I 

Then the solution of the constraints 

4.^, ie(I)+l^ a(N\^) (3.3a5) 

considered as equalities, is an extreme point of 1(1). 
proof j If r2(a(N\s)) = , 

then r^CaCuN^S)) = n-s. 

Hence, there exists a subcollection 

P(n\s) ,^a(H\s), i§(NV))j = n~s 

such that r^CPCnX,^)) = n-s. 
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Prom the earlier results, when the constraints 

-|- 6 ( I ) -i-0(l)+l _ \ 

, T 1 6 p(n\S) are made equalities, the v;hole set 

of constraints is equivalent to (3.3.3), (3.3.11) ana tne 

following set of inequalities. 


'RH(e(I)) = RH(e(l)+l) 

for all I with O < 9(1) < s-1 

. ^ ^i ~ ^ . ^ij S. / SH(o), for all I with 9(1) < 0 

i G I i ^ 3 

G I Rli(s), for all i v;ith 9(1) > s-li 

L 

Hence the system of equations , obtained by malcing all 
the constraints of (3.3,15), equalities, is equivalent to 
(3.3,3), (3.3.11) and 

2 X. - Z = RI1(0(I)) = RI-I(0(I)+1) for all I G a(l\s) 
i G I ^ i < D 

i,j G I ■ (5.3.16) 

Since rg(a(E\s)) = ^ 


(3.3.16) has at most one solution. show tho.t (3.3,16) has 
exactly one solution which is 

(2h+l) - a? 2a,- a- 

^ 4 . JL J. A 

~ ■ h(h+l) V ^i3 “ h(h+l) ’ 

for all 1< i, i,3 G(l\s). (3,3.17) 

It is not difficult to check that (3.3.17) satisfies (3,3,16). 
Hence the system of equations obtained from (3.3,15) has a 
unique solution, gi-^n by 


(3.3.3), (3,3.11) and (3*3.17) 
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The uiiiq_ue point so obtained, is an extreme point if and only 
if it is a feasible solution of it is merely a routine 

work to check that the point X gi-ven by (3.3.3), (3.3.11), 
(3.3,17) is the point (a,h) in Y* where 

s-I “ i G S and a^, i ^ S are predetermined 

and hence from theorem 2.1, X is feasible. 

This completes the proof, 

jL * easy to see that the restriction on 9(1) in 

the assumption (iii) of theorem 3,3 is to ascertain that both 
the blocks existing in E(S), for all 

I G a(N\S). Since relation (3,3,16) requires only the existence 
of the dominant constraint, it is sufficient to have atleast one 
of these two blocks where the dominant constraint appears. 

Thus it is sufficient to have the restriction 

-1 < e(i) < s. 

But from lemma 3,2 in order that the set of constraints are 
reduced to a set of dominant constraints it is necessary that 
for some subsets I, the consecutive blocks are 

equalities. Hence the hypothesis (iii) of theorem 3,3 can be 
relaxed a little as follows; 

(iii)' -1 < 0(1} < s ¥lea(N\^) 

and there is a subcollection §(n\^) oa(J^\p) such that 
|P(h\^)| = r^(p(H\s)) = n-s and 0 < e(I) < s~l, ¥ I 6 §(1^). 



2 : iletaods to generate tne e^treiae point by using 

ca.oxem 3,3 cojo oe summarized as follom/s; first a subset S mid 

integers ii > 1 o ^ n .\r 

/. a^, 1 e i'i\o are caosen. If tnese integers 

Sc-ois-/ use nypofcueses of tiieorem 3.3, tneii tue point 




.|.o, . . . ,a^, S V is an 


e..teeuie point of I'(I), Thus for rai possible c.ioices of 3 such 
oaat |sj > 2 and for all possible choices of integer h > 1, a. 

isiyiii^ tne conditions (i), (ii)? (i^.i) of taeoreru 3,3, a 
class of extreme points cmi be gener.-.ted mid it cmi be snown that 
such exl-reifie points ai'e -.aeiabers of Y* . 'Im.e class of sucn ex’treiae 
points willbe denoted by Eat (l}.Proui t leoremi 5.3* follm/ing 
results can be directly deri’ved. 

p;^or_ollaa-jr pu : 


For n - 1^2, 3, 4 all the extre-ne points of P(»') are of first 


ind • 


i I'roiii the examples discussed in section 3.2 the proof 
can oe directly derived. 


Oor pllar^^^ 3,2 • 


7p y y^OEoct (1) 


» Le S — 'lid — 0 lor all i 0 Tiien define 

Q:(i.i\S) = {{i>, {1,3 > J i < i, i, 3 6 Clearly for lx h < r-2, 

the hypotues--s of tneorem 3.3 are satisfied and' from (1.5,2; such 


a point is a member of Yq. 

- 



Similarly by talcirg S = = h\s, = -1. i e and 

defining a(H\s) = { { 1 1, 1 e > , it Is not hand to cheek that 
the hypotheses of theorem 3.3 are satisfied and from (1,5,3) 
sucli a point is a member of V,r. 

Thus so far a subset of Y* is classified such that the 
points are extreme points of the convex polyhedcon The 

cLuestion that still remains unanswered is that whether there is 
any other extreme point in V*. fhe foUaring theorem gives the 
necessary and sufficient condition for a point of 7* to be an 
extreme point of 

ghe o rem 5.4 : 

A point (a,h) G 7* is an extreme point of the convex 
polyhedron P(li) if and only if there existsa collection a(N) 
such that, 

( i) r2(a(N)) = , 

(ii) 0(J) = 0 or, -1 for all J e a(N). 

• Substituting the va,lues of x^ and x^^ in terms of a^jh 
in any of the constraints (2.1.2) we have, from 

1 ' 1 

\ h - E a. j ) h+1 - 2 a.f >0. 

1, 1 e J ^] ■(_ l e J t 

Por this constraint to be equality at the point J = (a,h) G 7* 
it is necessary and sufficient that 
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Thus X e (a,h) e T* is an extreme point iff there are SlS+il 
linearly indepenient constraints Bhioh are equalities at I. 
Since a(H) provides requisite numter of linearly independent 
constraints and condition on e(j), j g a(ii) ensures that these 
independent constraints are equalities, the rest of the proof 
immediately follows. 


The difference between theorem 3.3 and theorem 3.4 is 
apparent. The later is the generalization of the former in the 
sense that, in the former a subset S of indices is distinguisted 
for which the integers a. are fixed at the level 1 and on the 
0 tiler hand there is no such restriction in theorem 3.4. Another 
remarkable aifference between these tvo theorems is that there 
are examples to establish the result of theorem 3.3 whereas 
no extreme point can be found which satisfies the condition of 
theorem’ 3.4 but not of theorem 3,3. 

Eo cample 3.5 ; 

let N = { 1, 2,3,4, 5,6, 7,8 > 

S = {1,2,3,4 > 

Then the associated H(S) can be written 
_J 5 6 7 8 5 6 57 58 6 7 68 78 56 7 568 5 78 6 78 56 78 

0 


3 

4 
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XlBij clp- — a- ~l.a=:pQ — -I t-, V 

5 6 > ag 1 and n = 2. l<‘or any I C^'s\s, 

the blocks corresponding to -t 4---P- n 

b X ana i ^ are identified as 

the shaded blocks in fig, 3.4, 

Ihe condition (i) of theorem 3,3 is directly satisfied 

as 0 < h < s-2. from the figure the collection a(l\s) can be 
derived as 

a(li\s) = {{i>, {i,3|, for all i < 3, i,j g i^S}. 
le. a(N\s) = «5>,{6>,{7>,{8>,{5,6},{5,7>,{5,8},{6,7>,{6,8>,{7,8>} 

r 2 (a(N\s)) = lo, r^(a(l\s)) = 4. 

shall show that the condition (iii)' of theorem 3.3 is 
satisfies for such a(l\^S), 

It is not difficult to check that 

1 0(1) < s for all I e a(M'\^S). 

It is to be noted here that had the above condition been Tiolated, 
there would not have any shaded region for the corresponding 
subsets I . 

Y/e define 

P(n\s) ={{ 5>,{6},{7V,{8>}, thenr3_(p(hy)) =4. 

It is also seen that 

0 < 0(1) < s~l, for all I G P(k\S) 

Such condition ensures that for every I G p(ii'ys) there is a 

■'!p f' -I % 

pair of consecutive sha(fe d , ‘feb®, figure 3.4. 
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Thus the point 

(1, 1, 1, 1, 1, 1, 2 , h = 2) e V* is an 
extreme point of F(]J), 

It is also easy to check tha,t this extreme poiiit is not 
a memher of any of the 8 classes of ■'vertex points generated by 
Kounias and Marin *] . 

observed that in theorem 3»3j the number of 
constraints v/hich are made equalities in fact are more than 
the required number . It is not difficult to find 

exactly constraints from (3,3.15), which when made equali- 

ties, give rise to the extreme point, Tnese coiastraints can be 
picked up for the above example as follows; 

{1,2>, {1,3>, {1,4:>,{2,3>,{2,4>,{3,4> , {1,2,3 > 

{1,2,4>, {1,3,4>,{2,3,41, {1,5>,{2,5>,{3,5>,{4,5}, 
{1,2,5>,{1,6>,{2,6>,{3,6>,{4,6>,{1,2,6>, {7>,{1,7>, 
t2,7>, {3, 7>, {4, 7>, {1,2,3, 8>, {1,2,4, 8>, (1,3,4, 8>,{ 2,3,4, 8>, 
{1,2,3,4,8>, {5,6>,{5,7>, {1,2,5,8>, {6, 7>, {1,2,6,8> , {1, 7,8 >. 


It is merely a routine work to see that the coefficient 
matrix, with respect to the above set of constraints, is ncn- 
singular and that it is a subset of the set of constraints (3.3.15). 

It is to be noted that such choice of independent 

constraints from (3*3»15) is not unique and the choice of a(lf\S) 
is also not unique*' , Thus the're- ^‘''several ways to pick up 

the same extreme 
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point. Such situation refPT-o 4 -^ 

degeneracy. The high rate of 

degeneracy of the extreme ■. „ 

^ points malces it difficult to solve 

hP(N) hy simplex method. 

3.4 Ext (2) : extreme points of aeaond kind 


All the extreme points of ^(I) are not covered by the 
class Bxt(l). It is observed that there exists another class of 
extreme pomts which are not members of y^. me main aim of 
this section is to show the existence of a class of extreme ' 
pomts of which are not members of but can be generated 

from certain nonextreine point members of 7*. Assume that 
I = (a,h) e V* is not an extreme point of I'(h). let a(E) be 
the class of subsets of U corresponding to which the associated 
constraints are eq_ualities at X. Clearly from theorem 5,4, we 
have , . 




The purposs of this section is to choose some more constraints 
which are not in a(I) but which together with a(N) constitute 
the class a(]l) such that 

rg(S{H)) = 

and when all the constraints jn a(iT) are made equalities they 
give rise to an extreme poinit of P(N), 

For a given S ON and for given integers h, a-, i e IfXs 
where 1 <. h < s-1, let the subsets of be partitioned as; 
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71 = {icn\s 1 ~i < e(i) < 3 ^ ^ 
Ti^ ^ {icz it\s I e(i) < 

Tig = {I ci\s j e(i) > s>, 


where e(I) = h - z a- , 

i e I ^ 

From theorem 3,3, it is clear that if there exist 



such that 


rg(a(N\S)) 


and there exists a P(H\S)C a(l\^), jp(N\s)j = r^(p(H\s)) = (n-s) 

such that e < e(I) <^s-l ¥ I e P(n\s) then the point 

(1, 1, , . . 1, ••• 6 V* is an extrene point of F(N), 



but there exists a 


P(i\s) 7C such that lP(l\s)| = rj(p(Il\s» = n-s. 

Following arguments of section 3.31 - 3,34 it can be shoim that 
the whole set of constraints are equivalent to 




1 < 3, ij3 e S 




z 

i e I 



i RS(o) ¥ I e 71^ (3.4.2) 

= jaHj(e(i)+i),¥Ee tc (3.4,3) 
t ‘ tie % ' ( 3 .4 .4 ) 


’ "'f/ ; ''-''v’.:-! ';4‘ 
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By direct suBstitation of values of x-, . in terms of 

i''" Be shown that the point 

(1, 1, ..., 1, ...,aj^,h) e Y* 

satisfies the system (3.4.1) - (3.4.4) with (3.4.3) being 
equalities that is, 


1 X- 

i e I 


. 2 _ X. . = iiH(e(i)) = Iiri(e(i)+1), 

1 < 3 


ij ;) e I 


(3.4.3) 


for all I e 71 . 


The system (3.4.1), (3.4.2), (3,4.3)' and (3.4»4) defines 
a convex polyhedron which is a subset of I'(N) and hence any 
point of this convex polyhedron is also a feasible point of F(I). 
Based on this observation we prove theorem 3,5, 

Theo rem _^. 5 : 

Any extreme point of the convex polyhedron defined toy 
the system (3.4,1), (3.4.2), (3,4.3)' and (3,4.4) is also an 
extreme point of F(N). 


Proof : Obviously any extreme point of such polyhedron satisfies 

the equations (3.4,1) and (3.4.3)* and there are at least 't' 
number of independent constraints among (3.4,2) and (3.4.4) 
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show that this point is an extreme point of S(i;) it is neoessar; 
to show that the point satisfies' linearly independent 
oonstraints of (2.1.2) as equalities. Jollovvina the method 
adopted to prove theorem 3.2, it can be scorn tsat the set of 
constraints 



There exists at least one extreme poiat of (3,4.1), 
(3.4.2), (3.4.3)' and (3.4.4). 

Proof : Consider the comrex polyhedcon defined by tte 

following system; 

the system (3.4.1) and (3 ,4,3)* together with 

S Xr - 1 X. . < EII(e(I)) (3.4,5) 

i G I i < 3 

i, 1 e I 

f or all I 6 \J 71:2 • 

The p 0 

all the 
extreme point 




It IS Cleax that the system of inequalities (3.4.5) 

defers from the system ('3.4 P'1 a a \ -u a,, • , , 

(3.4.4) by the right hand side 

constants. 


It can be easily deri^d from a corollary of well known 
Resolution Theorem C 39 » page II9 ^ [_‘48^, CSO^l, 

L 31 J that any convex polyhedron which is defined by the 
system (3.4,5) by changing right xiand side vector, has at least 
one extreme point as long as the polyhedron has at least one 
feasible point. 


Since the system (3.4,1), (3.4.2), (3.4.3)» and (3.4.4) 
has a feasible point, namely J, it has at least one extreme point, 

This completes the proof. 

As the theorem 3.4 gives the necessary as well as 
sufficient condition for a point in T* to be an extren^ point of 
I'(N). Clearly any extreme point obtained by theorem 3.6, which 
is shown to be an extreme point of I'(R) in theorem 3.5, cannot 
be a member of Y*. Such extreme points are said to be extreme 
point of 2nd kind, let the class of extreme points of 2nd kind 
be denoted by Ext (2), , ’ ^ ^ ^ 

Example 3,6 : 
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Clearly rg(,t) = 

The reducedset of constrattits is given by 




Clearly abovB system has an extreme point as x 
H''nce the extreme point of I'(N) is 


~ I > ^ < i < 6 , x^ = 1, Xg = 1 

^ij “ 3 > 1 1 < D < 6, x^^ = x^ = I 



Corresponding X e V* differs from aboTe by x,j,q = 4/3. 

Thus we see that for an (a,h) e V* with = 1, i e S 
for some subset S if there exisisa collection P(h\^) such that 

r^(P(N\^)) = n~s and 

0 < e(I) < s-1 for all I e §(n\s) 

then either the point (a,h) e Y* is aii extreme point or, an 
extreme point can be obtained from (a, h) e Y^ which satisfies 

l6(I)^ l9(I)+l^ as equalities, ¥ I G She following 

theorem gi-ves the generalisation of this result. 



Theorem 3.7 
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gi-ren any (a.h) b T» which satisfies at least n constraints, 
which have rank of order one as n, as eciualitles then there exists 
an extreme point 2 of J(II) such that (a,h) and I have at least n 
common equality constraints. 


* J'or any extreme point X G i'(N)j let a(lj) be the set of 
constraints that are equalities at X. Hence 

r2(a(H)) = , 

r^(a(N)) = n. 

Let p(ir) C.a(N), jp(H)| = n and r^(p(l')) = n, Vfc claim that 
there always exists a set of (n+1) integers a^ , 3 ^, . . .,a ,h such 
that they satisfy the following 


I an = h or, h+1, ¥ I G P(I) 

i G I 


(3.4.6) 


Because the system (3 .4, 6) is a system of n linear equations with 
n+1 variables, by suitable choice of h, n integers can be 
determined. Hence the' point (a,h) G V* defined by these integers 
satisfies all the constraints J, J G P(N) for e(J) = 


h - 


a,. 


i G J 


0 J • 


Con-verse ly, let a(H) be the set of constrain-ts which ane equalities 
at (a,h) G Y*. If rgCaCl)) = , then the theorem boils 

do-wn to result obtained in theorem 3.4. let r2(a(h)) < ‘ '^3 ' » 

then by malcing use of th^efi/ii% -aorae more constraints can be 
chosen which together I". ^^aastitute a(I) sUch that 

"'iSf a(I ) are ' made 




r^(a(N)) 
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s qualii; ie s "they giy® 

y rise to an erfcreme point o_’ ]?(m). 

This completes the proof, 

2.5 ehai^^teriaatim of extrme joints of F<m 

In this section we will decior’Vio • j. 

scrj.be soaB ar tempts made 

0 and ail the extra* potote of the convex px-yhedron P(M). 

xe helleved that Ext (1) together with Ext (. 1 ) exhaust all 

the extreme points of I’d). rnii on ■ 

U^;. ihe following discussion is aimed 

ngthen the belief that all the extreme points of I’(n) 

can be generated from the extreme points of P(s), except 

possibly wilh one oxcoplion# 

^ ° “y of H having n-t element 

without any loss of generality we shall assume that S-t = 


f 1 9 4 ^ \ 11 ) 

) ,»».,n— , Clearly i'(]j-t) is contained in R ^ 

.(n--t)(n-t+l] 

Any point e H ^ can be defined similarly as in 

n(n+l) 

the case of Z e R ^ . 

i or the sake of uniformity we shall denote I e R ^ by z“. 
Clearly and haye some coordinates in common, namely 


the case of Z e R 


v^^-t n _ n-t p . , . 

» ^ij “ ^ij » 1 e II 


II -t. 


^ In this sectiaa it is sh(Mn that given an extreme point 
^ of I’(R-t) several exti©^ points of F(S) can be obtained 
such that these extrpme_|p|«ts, ^ye^ same values for tte coordinates 

. ' 7 " lA.' '*7 ’ 
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theory for generating extreme points of P(N), is of interest of 
its own. In [l3l^ tne class of points given in (1.5.1), 

IS obtained from the vertex points of 3?(H-t) by assigning the 
value zero to the new variables, arising due to the increase 
in the dimension. Here it is sliovm that such extension from 
I'(N-t) to F(H) can be generalised by putting several values, 
besides zero, to the new variables, lie consider tv/o different 
cases, namely when is an extreme point of first kind and 

when is an extreme point of second kind. In 3.53, we show 

that there is one extreme point of I'(Ii) which cannot be obtained 
from any extreme point of I'(I-t). ye believe that this is the 
only exception possible, ye give a conjecture in this regard. 
Such a belief is based on both theoretical and computational 
experiences. All the extreme points generated by random 
computation are members of either Ext (1) or, Ext (2), 

3 »51 : In this section we consider the possible extension 

of to extreme points of E(F) where is an extreme 

point 'Of first kind of E(]i-1). 

Since X^i^ is an extreme point of first kind, there 
exists a subset ScH-1 and integers h > 1, a^, i e such 

that satisfies the following set of constraints as 

e qualities, 

(j,^, ((,^+1, iQ(^), , for all I G a(Ii-l\s). 

Since S c. h- 1 ==> s CN, constraints of liP(N) can be 

arranged as R(S). Slk': vis partitioned into two 

■ bvb ■ i*-? 1 'i' ■ .. :: ' ' - 
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subsets such that the first one corresponds to subsets of N-1 
and the other one corresponds to subsets of H which invol^re 
the index n. With the help of this partition of columns of E(S), 
the array R(S) can be partitioned- to tv/o arrays namely, R' (S) 
and R*(S). Tnus R’(S) contains all the constraints of 
and R”(S) contains all tne constraints which are added to R’(S), 
to get all the constraints of II>(R) . 


Let us consider example 3.5, 


N = { 1,2, 3, 4, 5,6, 7,8 > , I-l = { 1,2, 3, 4, 5,6, 7 > , 
S = {1,2,3,4>. 

R(S) is given in fig. 3.4, 


The columns can be partitioned as follows; 
Columns of R’(S) are 


4- , { 5> , {6> , {7> , {5,6> , {5,71 , {6,7> ,{5,6,7>, and 
Columns of R"(S) are 

{8,1 {5,81, {6,81, (7, 81, {5,6, 81, {5, 7, 81, {6, 7,81 ,{5,6,7,81. 

We define any block of B’(S) as 


l'^’(S) = { UI 1 ®^C6j^(S)l 
for any I 

and any block of / 

R"(s) =Ml® 


( 3 . 5 . 1 ) 



( 3 , 5 . 2 ) 
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Ifence it is easy to see that there Is a one -one 
correspondence heteeen E'(s) and Ii"(s) because of the 
correspondence between (3.5.1) and (3.5.2). Ihus R-'(S) can be 
viewed as a copy of E’(S). Ihis partition helps to extend 
any extreme point of F(HJ) to an erfrene point x" of F(N). 

If x" is an extreme point of first hind of !(H-1) then 
e (As defined in sec. 2.3) 

^ ~ • * 'j^n-ii^) such that = 1, i 6 S. 

we assume that = (1,1 l,a^_,^^, . . .,ajj_ 3 , ,h) e 

The purpose is to obtain such extreme points of ^(N) so that 
the common components of and have the same values. 

That is, to determine the possible values of a^ such that 

(1, 1, ag^ 2 ,...,a^^ 3 _, a^,h) e V* 


either is an extreme point of P(I) or, gives rise to an extreme 
point of second kind of ^(U), 

is an extreme point of F(N~1). Ifenoe there are 

linearly independent constraints which are in R»(S). If 
is an extreme point of F(I) obtained by extending then 

at least n = (Hln-tli _ 

linearly independent constraints 
in R’'(S) have to he made equalities. Theorem 5.9 gives the 
condition for existence, of. spch n constraints in R*’(S). 




KI)+1 


substracting (3.5,4) 


If for any there exists at least one I r (i.!-l)\s 
such that 0 < KD < B-1 where ,(l) = e(I) -- then there is 
an extreme point of ^(h) such that 

_ Ji-l n n-1 
Xf - xt , X. . = X. . , for 1 < i < j 

P_roof i let for some I C (B-l)\s, o < *(1) < s-1 the system 
of equations obtained by making l*0)"(sj^ jHI)+l"(g^ equalities 


Xi + 1 Xj_+: 

♦ (I) ^ ® ^ 


2 X • ' " 2 X . =1 

i P m m . p T in 

1 e 1 e I 


(3.6,4) 




^i1 " ^ 

■>-& i e 1 

<i>(i)+i jei 


<f(i)+i 


(3.5,5) 


^(i)+i = 5:^1) O at forke s\i^(y. fhen substituting 
the value of in in,.., (3. 6.4) - (3.5.5) and 
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Simplifying, we get x. 

kn ir(h:+r) (3.5.7) 

Since (3.5.4), (3.5.5) hold true for all 1 ■ p . .qv . 

aa.± Q 6^^j^(S) and 

^(!>(I)+1 6 ^(j)(l)+l(S), nence (3.5.7) holds for all k e S. 

i!ow using similar argument as theorem 3.2, substituting 
the values of (3.5.7) the whole set of constraints in a column 

of R"(S) can be reduced to either a single constraint or, a pair 
of constraints as follows. 

for any column I of Il’‘(S), 


i e I 


That is, 




1 


X 


i e I 


in 


< RH <J(I) = RH(4i)+ 1), if 0 < ^(I) < B 
1 I^(0) , if ^(I) < 0 

1 ^(s) , if (j(l) > s. 


•— ' n(h+'l) ' ^ kfO_£^(I)_<s 

ay^(a^+l) (3,5.8) 

^ “T(1mTJ » *^(1) < 0 


a^(a^+l) + 2a^ 


(3.5.9) 

, if <>(I) > 3 (3.5.10) 


- “~hTIi+r) 

Using the results of theorem 3.3 and 3.5 we conclude the 
following I ' , 

( i) If there exist (h-^) number of subsets of (h-lj'^S such 

that for all these ,,subs#ts iiptber; 4(1) satisfies 0 < 4(1) < s, 
and, when the set of inegu^ities' (3*5.6), associated with these 


n-s subsets, are 


anique solution, then ' 

with values 


rise to a 



the -first kind 
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_ ^2h+l) 2^1 % 

= Tj;m?l 1 e (E-1) 

and. rest of the coordinates are same as 

( ii) Otherwise j the situation that at least for one subset 
Ij I G (U— 1)\^, 0 £, s— 1 ensures that by selecting sohb 

more constraints from (3.5.9) and (3.5.1o) we can get (n-s) 
subsets of (E'~l)\s and ¥vhen the corresponding constra.ints in 
(3.5.8) ~ (3.5.I0) are made equalities tiiey give rise to an 
extreme point of the se cond kind. 

In fact the above two conclusions can be directly 
derived from theorem 3.3 and theorem 3.5. 

Hence the proof is ccmiplete, 

5 .52 ; Here we discuss the possible extension of any 

extreme point of second kind of HCH-l) to an extreiB point 

of F(H), Let (1,1, . ^ point in PCH-l) 

from which we generate the point G Sxt (2) as in section 

3.4, Hence this point satisfies the conditions given in 
theorem 3.6. The aim is to find the values of a^ such that 
the point (l,l,...,l,ag^;L»***»^-.l’^’^^ ® V* either is an 
extreme point of I'(N) or, generates an extreme point of second 
kind of P(H) . The common coordinates of which have 

same values are given as^.., : ■ 



5 



where e ie defied in theorem 3.5. In view of theorem 3.3 end 
tneorem 3.5, liie theorem 3.9, conditions a in terms h, a., 
i e N-1\S can be derived* 

Thus such a point can be viened as the extension of 
the extreme point X® of l'(S) to P(l), where x® is 

^i = BTT ’ ^i;i = h(l+rj » 1 < i < 0 < s. 

Clearly X® is an extreme point of first kind of P(S), 


^*.55. * have shown tnat there are some extreme 

points of P(N) which are obtained as a result of extensions of 
Ext (1) of P(N-t), t > 0. Here we give an exceptional case, 

We give an example of one extreme point of P(N) which cannot 
be obtained as a result of such extensions. 

let X^ G P(N) be defined as 


n = sir > nj = . 1 < i < a < n. 

Clearly G Ext (1) and x'^ G W*. 

In order that x"^ is an extension of any extreme point 
of P(S), S it is necessary that the point X® defined as 


x| = x^, x®j = x^j, i < 3» fj3 e S 
is an extreme point of 1(3} ^ 


. X'*'/ " • 

^ -‘I ' . 


But from theorem 3»1 




prsl for n <'h < s-£ ahd hence Z® is not an extreme point of 


of l'(S) for 0 <'hk s-e 
PfS'). 




:^*isa«saiia 





We elleye that except for the ahow particular case all 
the extreme points of eia he generated hy the method 

given in 3.51 and 3.52 from the extreme points of F(N-t), t > 

and all these observations are summed up in the following 
conjecture. 

Conjecture 3,1 


Per any extreme point of F(i) there exists 

subset S such that 2 < s < n and 


s-T ’ ^ij -[T-TJTCT)’ ’ ^ ^ & S 

3,6 Extreme -points in 7^, 1 ^i ^8. 


fo strengthen our belief we show that the classes of 
points . 7g, Vg, Y^, Vg, which are not in Ext (1), Ext (2), 

are not extreme points of i’(N) 

theorem 3 .In 

The vertex points Y^, Yg, Yg, Yr^, Yg defined in 
( 1.5, 4)-( 1,5.8) are not extreme points of the convex 
polyhedron of E(N). 

PJ^oof : In theorem 2,1 part II it is proved that the members 

P' 5 > P'gj <^3n be, expressed as the convex combinaticin 
of Y* , So it only reitei^'.it'a .show , that nene of the points 
of Yg is extreme p from definiticai 1,1, any 

to some linear system of 


Hence eny point which le solution to tnla system can not 
have irrational coordinate values. iny point i e Vo vvhich 


is defined, by (1,5 . 8 ) from the set of numbers b^jb,., 


. b 


out of which at least one 
we haAre 


n 


b. 


is Irrational then from (1.5,8) 


at least one coordinate of X as irrational value , and hence 
such X can not be an extreme point of the system ^(N) . 

Ihus it is sufficient to check only those X e V, which are 

d 

defined by the set of rational numbers b., 1 < i < n, 

1 ***** 


iny such numbers are of the form 



where are relatively prime, q_^ > o and h defined 

as h = ICSi J . 


¥/e shall show that when h ^ 1, the associated X G Vg can 
be expressed as the convex combination of the points of ?* 
and if h =' 1, then the point X 6 Vg is in the con^rex 
combination of two feasible points of F(H). 

Case i If h 7 ^ 1, we determine two points of Y* as follows 


l 03 


( 2 h~l)a.-af 

-V I f — 1 1 


TTfErrj- 


5 x! 
ID 


2 a^ 8 ^ 


Now X is in as follows 

h +1 


Y — f _i. • 

- TTCT' a' + iL' 


To see this 


2 h+l , h -1 

'i * 


•ar + fr =^1 


2 h a^-af 
72' 


? 

® E- - 5^5 = Sb.- b. = 

2 a_. a 


-I iH"l I h— 1 , , 

2¥” ""ij + '2Tr ""ij = 


h 


'H = 2 “i“f "^ij 


Thus f or Xe V,Q, there are two point 1 ', l’ ' e such that 
X e g[;x',x -Hence such X can not he extreme point of 
P(N) . 

Case ii If h = 1 i.e., if all h^ 1 < i < n are integers 
for the point X 8 Vg under consideration, in this case 
case X' can not be defined, for such cases it is shown 
that for any such X 8 Vg there exists a point X' 8 V* s ’ch 
that the point 

Y = X‘+ x(X-X'), for any X > 0 
is a feasible point of fCNl-#, Hence IS Since 

' ' ' : ' -y ’ --4 , '• ■ 4 - 

h = 1 ==> b. = a-, 1 < i", f;,h» X* 

^ ^ :-^:r 

■ 1 :' 


t# 


:x:. 






m^m 
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1 ov/ 


b, 


Yi - + x(x.-xp = 

^ir X(0ib.) 


^ 3 ' 


for any arbitrary 

constraint 

J CH 

2 


yf 

- 2 ■ 

7 ^ - 

i e 

J 

X 

i < j 





i, jCJ 


= 2 



2 x‘. .+ 

X{ 2 

i£J 


X 

< j ^ 

iej 



i. 

fDSJ 


= 2 

X' 


2 x'. . 

+ s { 

ieJ 

X 


i<j 

2 




ij jeJ 



b.-b? 

i 1 


ieJ 


- £ b . b . } 

i< 0 3 

i, joJ 


ieJ 


Since for any integers b^, 2 b^- {£ < o. Hence, 

for X > 0 

2 y.- £ y . < 2 x! - 2 x! - < 1 

ieJ i<3 ieJ ^ i<3 

ij jeJ i, 3SJ 


for 

X‘ 

is 

a feasible point 

==> 

Y 

is 

a fea,sible point 

==> 

X 

is 

the convex combination of tv/o feasible point 

==> 

X e 

^8 

is not an extreme point. 


This completes the proof . 

3. 7 Adgaaenay struatufe of extreme points 

In the section,,, 


aspects of the adjacency, „ 

. , - V % •- 



is ;»ade to study seme 
the e:rbreme ootnts 
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‘2vi0 extreme points x', X" of i(N) are said to oe 
adjacent if there exist (N) C y' (H) ? cc'J(i) C y”(iJ) 
where v’(h), y''(N) are the set of constraintswhich are 

equalities at X',X" respectively, satisfying the following 
relations . 


'lie discuss here the adjacency among a pair of distinct 
extreme points in Ext (1) satisfying the following conditions 


he a subset S, h > s-2 such that 


where 6 . i 


¥e here assume that there exist 


such that a!’ 


first kind, it follows that 
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~ Jt ""I i e cc'(iAs) 


-1 < - _ E aP< s ¥ I e a'«(iAs) 

IGI ^ \ 


(ii) There exists P ‘ (l\s) C a' (A^) j P' ’(M\3) C a’ »(jjy ) 
such that r^(p»(N\s)) = j{3>(h\s)| 

= r^(p"(h\S)) = jpti(H\s)j = (n-s) 


= r. 


0 < h - a| < S~1 ¥ I e p»{h\s) 

X0 1 

0 < h - aj_' < s-1 ¥ I e P'»(h\s). 

10 I 

how let us define the following collections of subsets 
of h\s . 


= 

{I 

c 

1\S j 

1 2 

a! 

- 2 

a! 

’1 = 

0 } 





iei 

1 

,fel 

1 


r: 

{I 

< 

h\s J I 

1 ^ 

a| 


a!' 

' j < 

1> 




\ 

iei 

1 

iei 

1 



r: 

{I 


h\s ; i 

2 

al 

- 2 

ar 

'1 > 

1) 





iei 

J- 

iei 

X 



y» 

"o 

a 

nd 

Yi 


■ = 

h\S , 



In order to find the condition that 


YI v! I 


are adjacent in 


^(h) we first prove lemma 3 .6 which is useful for the purpose 



Lemma 3 .6 


show that 
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\{e know that ' (N^S) ) = n-s 

n P'(N\S)) = n-s 



Hence 


(3.7.1) 


iei ■ iei "'i' ^ 

Prom (3.7.1) clearly the system of eq^aations(3.7.2) is 
equivalent to the following 


a = a^ * ¥ iei\s . 

Hence is a contradiction, for we have a,ssumed that at least 
for one ieHN^, aj_ aP . 

Hence P’(N\S)) n-s and the proof is complete. 

Th eorem 5.11 

If r^(YQr) p'(N'\^)) = n-s-1 (3.7.3) 

and cc*(il\S)) = (3.7.4) 


then the points X’ ,X’' are adjacent in P(M) . 

Proof ; Since ja‘(N\S)j = r2(a(N\s)) = , 



Again frcan (3,7.3) we 
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Clearly such I ^ hut sjuice l e p’(I\S} C a'(K\S) 

A. 

==> I e ^ , from (3,7.5) 

Hence | 2^ al - a'*. I = 1 

iei- ^ iel ^ 

==> le'(i) ~ e' '(I) 1 = 1 

where e’(t) = h- E, a! and 0"(I) = h - E, a!’. 

iSl iel ^ 

Let us assume that 0'(i)-e**(i)=l 

==> e''(i) = e'(i) - 1 (3.7.6) 

Irom sections 3.31-3.34, we kno?/ that the constraints 
I® I® are equalities at the point X* and 

similarly are equalities. Fran (3*7.6) 

we have q® ' * . 

Ihe set of ccsnstraints v/hich are equalities at X* are given 
hy, interms of blocks of E(S) . 

^h^ ^h+1^ I®'<^^), ¥ I e a«(I\S) (3*7.7) 

Similarly the set of constraints which are equalities at X’ ’ are 


<l>^, a«‘(il\3) (3.7.8) 

Using the relations (3.7 .3)-( 3. 7. 8) we can writs some 
constraints Y?hich are -equalities at both x* , jJ ' as 
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J- 

From the definition 3 ? i-i- 

It suffices to show that the above 

collection of subsets has rpni- t n(n+ l') . 1 

^2= -Y— ^ -1 and wnen 

either ie(Hl or, " (D ^3 I 3 

Eepeattog the step by step argument used m section 3.3 the 

above facts can be proved. lo avoid unnecessary repetition, 
we omit the rest of the proof. 


Ihe following exaraple illustrates the above result. 
Example 5 .7 

Let us take the same example as example 3,3. Let 
be defined as 

a,j,=2, ag= -1, h=2 (Fig. 1.4) 
and X’’ be defined as 

ag~lj ag Ij a^— Ij ag=Oj h =2 

a‘(ir\S) = {{i}, {i, 3 }j ¥ i < 3 , i ,3 e K\S} 

and a' »(N\S) = a»(X\S) 

p’(N\S) = {{5} , { 6 ) , {71 , {? 8 } } 

Yq = { { 5} , { 6 } , {78}} 

= 1\S' and Y> = $ . 

(p‘(X\S)) 

and 
Hence 



Ho 


Mce example 3.3, we can choose suitably ife+ll 
oonsWnts lor each x'.IM ceperately end show that they 
have -ift-L _i constraints as egualltles In common . Ihe 
set of constraints which will define x< is 

<12} , (13) , (14) , (23) , (24, _ (. 4 , ^ ^ ^ 

(134) , 1234) ,115) , 125) , , 35 , ^ ( 45 , ^ ^ ^ 

(26) , (36) , (46) , (126) , ' ( 7 ) , ( 17 ) , ( 27 ) , { 37 } , 

(47 ) , (178) , ( 278) , ( 378 ) , (478 ) , (1278) , (56) , ( 57 ) 
( 1258) , (67) , (1268) , (1238) . 

The set of constraints which will define x’ ' is 

(12) , {13} , {14} , { 23 } , { 24} , { 34 } , { 123 } , {124} , 
{134} , {234} , {15} {25} , {35} s { 45 } j { 125} , {16} , 
{£ 6 } , {36} , {46} ,{126} , {17} , {27} , { 37 } , { 47 } , 
<l^y f {178} , {278 } ^ {378} , {478} , {1278 } , {56} , 

{57 } , {1258 } , { 67 } , {1268 } , {1238} . 

V/e can see that except the constraint {7} rest of constraints 
which (iefine X’, also define X” together with the constraint 
{1,2,7} • Hence hy definition 3»3, 

X', X** are adjacent, 

Remaric ; In order to show the adjecency among any given pair of 
points suitable P'(l\S) a^d,, should be chosen, lor 

instance, if we take, ; 

X''=(aiH,Xa2=4:t%'^/\A'a5i=l, ag=l, a,^l, ag= -l,h= 2 ) . 
Then by choosing a different p*(S\S) it can be shovaa that X’ 


and X‘ f are 
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Chapter iv 


G^BS 

4» 1 Introduction 

This chapter deals ■with sorne special cases for which the 
best Bonferroni lower bound can be coraputed. 

Section 4,2 gives some useful results concernii’^ the 
special cases discussed in section 4.3. Section 4.4 deals 
with the condition under which a specific extreme point will 
be optiQial. 

4,2 Elementary results 

This section presents some results concer-ning the 
constraints that are equalities at the optimal point of If(N). 
Though all the results are not directly used for computing 
the bounds some of them are irportant in this respect, 

Lem~ma 4.J. ; 

If for some i, 

\ 

then Xj^ = 1 at the optimal 

Ir ppf ; \/e prove 

that X* be the 




.;nV: 


point of 




Let US ctHisider 
I X i < n, 



se show that there Is arother feasible potet sues that 

ip = 1 and the -value of the objectlTC fimotlon at 2 Is not 
worse than that at X*. 

Let d = 1-x^, 'then clearly d > O. he fine the point X 
as follows, 



(X* .3^ 

/ 

' - ( 

^3k » 

3 7^ i, k ^ i 

^3 = 


’ * 




(^x*+d, 3 = 

i 


i“lCj 3 7^ !• 


¥e claim that X 0 i'(N). let us consider any censtraint N 

such that i 0 J„, then 


Z 

k e Jj. 



z X. . 

k < 3 

k,3 e Jj, 


z 

k G J . 


^k- 


z 

k < 3 

k, 3 0 J, 


* 

^k3 


+ d -* (r-l)v 


(4.2.1) 


If r > 1 , then d - (r-l)d < 0 - 
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Stooe the oonetralhts which ao not the index i 

are not affected by the change made from I* to X, (4.2.2) _ 
(4.2,3) shov/ that the point X G I'(E), 

the difference in the value of the ohjective function 
at X and X* is given by 

P X - Pd* = dCp, - I s ) > 0 
Hence P^X > p-"X*. 


'Ihis completes the proof. 


lemma 4,2 : 


If for some subset S and i ^ S, there is an integer 
k such that 


k-1 


P - - Z P . . + 

1 3^i ^ 3 G S 


2 5ii>0 

G S 

then at least one of the constraints from 


(4.2.4) 


k-l 

U (6t(S)U {i}) 

t=0 

is equality at the optimal point, 

1 y po_f : let X* be the optimal point. Assume that none of the 


constraints of 


V/ Ul)- is - '^.quality at I*. Define 

tX) l' 


a point X such that for some 


r ' ' -vh^"; t* • ' . 

I ome 

i/'. ’ ^ dti: 







t=o 

In order that the point X to satisfy above constraints the 
value of d is restricted as follows, 

d = Min ! 1 « 2 rf + E x* - "1 

hh.J 


Min ^ 

L j e J 3 

j e W ( 6 ^(S) u<i>) 
t=o 


E 

h 

< Dg e J 

(4.2 ,6) 


lor any d > O satisfying (4.2.6) the value of the objective 
function increases from x* to X. Moreover (4.2.7) ensures that 
one of the constraints of 


lemma 4,3 : 


k-l 

(6-^(3) {i}) is equa,lity. 

t=0 

If for some i, p. - E P.. < 0, then the 

1 n-2 13 

constraint associated with is an equality at the optimal 
p oint , 


Proof : We follovi? the same method of proof as in lemma 4.2. 

let X* be optimal point and p. - E P - . < o. Define X 

1 n-.; 13 

1 

d, 3 / i for some d > O 
xt , 3 / i, k ^ i. 


as , 




A 

^i 

= xj - d, 

A. 

md = 

TC* 

^i3 

and 

^k "‘■k* 

A 

X = 

3k 

X* 

^3k 


It is easy to see that as d increases the value of the 
objective function increases and the point X remains feasible as 
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long as the constraint is satisfied. Hence the maximum 
increase in the ohoective function is possible by the maximum 
possible value of d, i.e. when t satisfies the constraint J 

n 

as equality. 

Ihis completes the proof. 

4, 3 Spetyial oases 

In this section we give some of the special cases for 

which the Bonferroni lower bound can te oalouleted by direct 
inspection. 

theorem 4.1 deals with one of the simplest cases which 
is ignored in earlier ¥/orks. 

Theorem 4.1 • 


is given by 


^i ~ ^ij ~ then the optimal point of Ii>(N) 


. ^i " ^ij = 1 < i < D < n 

and the exact lower bound is 
n 

I 1*1 - 2 p . . 

1 < i < 3 < n ^3 

! ^rom lemna 4.1 we know that = 1 for all i at the 
optimal point. Subatitutlng the values of x, in IPCl), „e have 

Max ,1 P. - B p, „ 

1 _< i < y < ij id 


2 

i < d 

i 5 3 G J . 


^ij<i 


1 < r < n. 


such that r - 
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2iae solution to the aho-ro linesu: program cmi be obtained 
directly as follows. Since P - • > o, 1 < i < i v >- ^ - 

> 1, the optimal point is given by 

Phis completes the proof. 

Theorem 4.2 gives a generalization of the result by 
Gallot and Kounias E29J. They show that for n = 2, 

n l"! = Pg = the lower bound q- j I 3 exact, ye give a 
wider class of points for which this bound is exact. 

Theorem 4.2 : 

It Pij = Min (P.,pp, 1 < 1 < 3 < „ raen P^ Q- p i^ the 
optimal value of liP(li), 

ItoM : IVe shall show that there Is a feasible soluticsi T° 

^ I£“(lf) such that the value of the dual objective function at 
Y -3 P Q- P. Since there fa a feasible point In »(») euoh 
that the objective function at thle point has the value P* q- p 
from the theory of duality we have p® q“ p in -, „ - ■ , 

of dual pair of linear programs. 

Pet Pi > Pg > ... > p^. 

Consider the subsets 
1 ■ 

^ ~ ~ = { 1 p 3 1 ^n . 

a'e define the dual variable Y° as 

y°(j‘) =Pi -P^^^ l<t<n-l 
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y°(J^) = and y^(J) =0 ^ J , 1 < a _< 

It is easy to check that 

£ y°(J) d(J) = P and y°(J) >0, 9 J C h 

J C N 

Hence y° is a feasible point of I3P(H). 


n 


Ihe value of the objective fimction at 1° is given by 


n 


Of -,k^ _ 

■£, y (J ) = Pi* 

k=l ^ 

It is also seen that q"! is given by 1,0,0, ... ,0 ]]] • Hence 

Phis completes the proof. 


P^Q-P =P2_. 


Remark : In order to ccxapute the best bound for such cases 
the following relation can be made use of 

P^Q"P = Max P.. 

1 < i < n ^ 

In theorem 4,3 v/e consider the case when all P^ are same 
and all the P. - are also same. Ibis case is aiso treated in 

ClsD# C29l3t CsiD , Here we derive the result 

using duality theory. 

Ih eorem 4,5 s 

If P^ = P^, P^^ = P^g, 1 < i < D < n, then the best 
lower bound obtained is 

E+I ^1 ktk+l) -^12 


(4.3.1) 
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( n—l ) P 

where (k-1) is the greatest integer less than — , 

^1 

P£oof : She proof is similar as that of theorem 4,2, It is 

shovm that the -value (4,3.1) is taken on by the primal and dual 
objective function at their respecti-ve feasible point. 

I>efine X e P(H) as, 

^ ^ij = kfkliy • 

Ihe dual feasible solution 1° is defined as 

y ( Jjj;) = y 

for all G 6j^(N) 

^k+l ® ^k+l^^^» where, 

r - - ~ 

"^1 - n~X““~- LHTf - 1*12 


"^2 - nflj 


k 


p k-l p 

^12 "iCT^l 


Clearly ~ ^2 2 0* It is easy to chedc that 


^k ® 


d(Jk) + rg 


‘^k+l ® ^k+1 


(I) 




Hence Y° is a feasible point of IP(H) 


Ihe value of the objective function at 1° is 


n 






n 




£-1 ^2 


(4.3.2) 


Simplifying (4,3,2) we have (4.3.1) 
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The following theorem is a generalisation of e^irlier theorem 
4.3. Though the result of theorem 4.3 is also obtained by 
Kounias and Marin such generalisation becomes possible 

only when it is proved by the help of duality theory'. 


Theorem 4.4 : 


If P^^ - P^g, P^ P^ji for all 1 < 3- < 3 < n. 

then the best lower bound is given by 

]ic+r^ klk+l) n k(k+ll ^In 




where t-1, k-1 are the greatest integers less than 

■^n ” -^In 


Proof : The dual feasible solution can be defined as, 

y y ~ ^ 2 ’ ^ ^ ^t+i^ ^ ~ ^3* 

y°(2^t+iU {31}) = r^, 

for all Tj^ G ® 

T^ e 6-(-(h-l), T^_^^ e 6^^^(H-1)5 

and, y°(J) = 0 for the remaining subsets H. The valu^ of 
r^f rg, r^, r^ are defined as 
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n-1 


n- 




C 


t~l 


'n 


n-1 


In 


] 


n-1 


^ 4 : 


t-1 


^in - ^ 


n 


It is easy to see that is a feasible basic solution of 
IJP(l'). It also provides an optimal point of IS’(ii), because the 

! 

value at this point is same as the value of the objective *, 

function of liP(N) at the point, | 

I 

‘ 

^i » ^ij ^ Ei[in7 9 ^ S. ^ ; 

- i ( t +l)(t+2) _ ^-t) 1 / w ^ 

^n ^ » ^in k(k+l) » 


Remark : Such cases as discussed in theorem 4,3, 4,4, can be \ 
extended to more general cases, but the simplicity v/ith which | 
the best lower bound is calculated is not retained. Hence, j 

further generalization will not lead to sinple computational i 
procedure and hence are not of much use. The possible extension! 
of theorem 4.4 can be thought as, ! 

f or some S C. IT, i 

^i "" ^1» ^ij "" ^12» ^it ^ ^in» ^ ^ i, j 6 S, t gf S j 

^t ~ ^n* ^tk ~ ^n-ln’ t < k, t,k ^ S. j 

! 

Ror such problems, the condition of optimality can be generalisei^ 
following theorem 4,4, Ihis problem is treated by Kounias £]30 3!| 
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4,4 Optimality of specific extreme points 

In this section we discuss conditions under which certain 
special kind of extreme points are optimal. In order to 
determine such condition, we malce use of the following results 
of duality theory. 

If X G ^(N) is an extreme point of I'(N), then there are 
independent -vectors d(J) such that d(J) X = 1, In 

m 

these vectors d(J) correspond to columns of A associated with 

the dual variables y(J). Let B be ihe square submatrix of A^ 

obtained by p idking up these colu-Qins of A"^ . If B is 

nonsingular then B constitutes a basis of I!P(N). Brom the 

well known theorem of duality theory, if > 0, then X is 

the optimal point of LP(]I) , To show that B is nonsingular and i 

® ^ > 0, we adopt the algorithm suggested by Murty [^39 2]* I 

Since the rows of B are coefficients of* components of X in LP(N),i 

we denote the rows of B as E. corresponding to component x- of X ^ 

1 1 ! 

and E^^ as corresponding to component x- - of X. [ 

Ihe matrix B is written in a tableau and in the course of [ 
the algorithm the rows of the tableau are altered such that at | 
every step each row in the tableau is a nonzero linear combination 
of the rov/s of the initial tableau. If at a particular stage 

all the entries of any row or column are zero, then the matrix i 

IS singular because Ihezero vector is shown as the linear | 

combination of the rows of B. Otherwise the algorithm stops I 

when by successive row operation the tableau contains the I 
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Identity matrix. If any elementary r„„ operation la represented 

algeloralcally am. E. + M . which means that E times of row E. 

is added to row E^, then the rows E,,E. . of the final tableL 

Obtained after successive row ooerat-i nnc 

upcra-cions, car oe represented. 

as linear algebraic function of R. r 

1* ij* 

let as first consider the extreme point x’ , where 




sSx ’ 


n 


( 4 *4 


21-i- ' In-llXnl^J ^ < i < j < 

Ihe constraints which are enualitiAQ r,+ vt 4-. 

equaxi-czes at X', are tae members of 

the family 

li dy, iCN), 

{4.4.2) 

whose oardtaallty is £Lh+ll. let B be the submatrix of A® 
obtained by choosing the column associated with 4„_g(10 o & 

The orem 4.5 : 

B is a basis of IE>(N). 

proof i Ihe aim is to show tnat 3 is ncnsingulor. Let 

®l3 1 r < 3 < n denote the rows of the tableau after BOme 

successive row operations. Sfe define £, E. . in terms of 

original rov/s as follows ^ 

21 


»i 


= [((h-3)Hi - .I,E,p H 


u 


n 


= sir V rn-^7fn^7 , , ^ ^ ^ 


^ t < 3 < r 

(4,4.3) 

Ka .-fi.-Xi.+B. . 


iP]M) 


(4.4.4) 
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It is not difficult to check that such R. , R. - transf ormi the 
matrix B to an identity matrix and hence 3 is nonsingular. 

This completes the proof. 

In the tableau discussed earlier, if the vector P is 
written as an additional column then the set of row operations, 
that converts B to an identity matrix, transforms the column 
associated with P to in the final tableau. 

I heorem 4. 6 j 

1* is optimal if and only if 


((n- 3 )P,-^E^PlP -5:3 ((n-3) ^ ^ ^ ^ <> 


(4,4 .5) 


nil ^ ^ t £ p _p > Q 

^ t=l ^ (n-l)(n-2) 1 < t < k < n ^ 3 iD *' 


(4 .4.6) 


i file R^, Rj_^, 1 < i < j < n defined in (4.4.3) - (4,4.4) 
transform;, the component P^ of P to left hand side of (4,4,5) 
and the component P^^ of P to left hand side of (4,4.6). Hence 
the condition (4,4,5) ~ (4.4.6) is same as > 0 v/hich is 

the sufficient ccndition for optimality. 


The ’necessary’ part can be derived from tne fact that 
bhe B obtained from X’ is uniq.ue for the point x’ and this 
LS not the case with arbitrary extreme point of P(H), 


Shis completes the proof 
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We next determine the sufficient condition that the 
extreme point x” is optimal where x” is defined as 

^i = eIt » ^ij = EP+rj ^ 1 i < 3 < n. 
for 1 < h < n-2, 

The constraints which are equalities afe the point X” are the 
nemhers of 6j^(U)u whose cardinality is 



^e choose the following subsets from S^i'SJSj 

nd form a basis of lE’(N). 

Ihe subsets are chosen as follo?/s 

{1,2, ...h} ij{j} ¥ 0 0 {h+1, ...,n} 

(1,2, . ..,h,h+l,h+2 > \ {i} ¥ i e (1,2, ...,h} 

■ d { 1,2, . . ,,h-2 ^ y{i, 0 } i < j, i,3 6 {h-l,h, . , .,n> 

{ 1,2, ,..,h-2> U ^h-l,h} U (k) V'b> 

¥ t G { 1,2, . . .,h--2 ) 
k 0 {h+1, . , .,n} 

{1,2, . . .,h-2>'\{t} }^{h~l,h+l} 

{l,2,...,h-2}\{t> ij{h,h+l> 

¥ t 0 {1,2, ..•,h-2} 

I {1,2, .. .,h-2>y {h-l,h,h+l,h+2>\{i,3 } 

■ i '"C „ J ^ i ijf J 0 C 1 j S j » ,i • j: ii»— ^ 1“ * 
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let us denote by a(i,), the family of above !l%li 

- - followin. theonem . determine o^dition! tbi: I,; V.,. 

^®asible basis of IS>(it). jhat ic, 4 - ^ 

, ^ ^ ^ if 3 is thp 

submatrix of A-^ chosen with respect tn p 1 

, respect to elements of a(I), then 

B IS nonsingular and > n. !i>hnc! 

sufflo-p 4 - ^ condition gives the 

sufficient condition that x" is ontimoi • 4 . 

optimal point of liP(I). ihe 

value of the objective function at z” is 
2 ^ 2 

^ i=l P,, = z* (SB.V\ 


< .1 < n 


^ 5, i < j < n 


^i “ ^ P • • 

j^i 13 

p !' = p ’ 1 S f 

•^1 ^i ^ 2 p’ 


=1 3 


theo rem 4,7 • 


^ i = 1,2, ...,n. 


= |a(N)j == aln+l}. 

and a(H) yieido a feasible dual basis if E . , p 1 < i < ■ , 

eatisfy the following set of inequalities/ ~ ' “ ” 


ri' > 0 


l,2j,,,jh 


< 0 . 

P . . - P . ~ p 
13 v^i 


r> 

L 

p!* 

> 

i=l 

1 


h 



2 

pt’ 

> 1 

i=l 

1 


i 

= h+2 

+ z' 

> 

0,: 


3 “ i< 
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n 


n 


D=1 


n 


- 2 p > 0, i = 1,2, . .,,h-2 

j=h+l ^ . 


® i ih-l-* Ih-^ ih+g-^i+f h+2+^ ' 

1 

ra ^b-lj + ^ j 2 Oj 3 = ii+3, . . .,n 


n 


^ > 0,i = 1,2, . . .,h-2 

;j=h+l 3 " » 


■i 


+ - ® * 3 “ ii+3, •. .,n 


h-lh 


1 < 1 1 


n 


£ P . . 

3 < n ^3 


^1“ 2 jp • - 

i— 1 h-1 _< i < 3 < n ^3 


ii-2 




st ’■j^ H . . . ^ 

^ 1<1<3<21 


> 0 


h-JL 


n 






n n 

2 P^^+ £P^ + £P- 

3 ' 


h-1 < i < 3 < n 3-3 Is "^23; 


> 0 


s = h-1, h, t = h+1, h+2 


^h+lh+2" [^__g < i^< 3 < n^i3 ^^ih-l+^ih) + ^h-l+^] 


n 




> 0 


P££of : The sequence of row operation, that transforms 

tableau containing B and. P to I and B ^ respsctively, ia 
given by 

q=(h-l) (e,- ^ .2, E,.- sir 1^ R.* R„ I 
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(^~ 1 } i — Rj ' ^ 


h 




^ e:t - 5- ^1 ~ 2 (Rt - J--r Z R ) 
0/^1 t=i ^ n -1 Rj -^) 

h ^ 

*'■ ETT 2 R. ~ rxr-^-m ^ 

5=1 5 n>+i)Xh-i 7 ^ ^ ^. . 

= ( b - l ) S - Ej ^ + D jj , / 

I ¥j ( 

^id “ Rii - fii - E ,.+- S - E 


^ hrf t=i t ^(^>+17 1 < t < fc < „ 


R . 


n 

ih + 1 ' ^ Rit -- 


n 


t=R+ 3 ^^'^ 


, n 

■ih+l^^i+Rh+l ~ R-; 

j=h+l 


+ 


2 ? 


^ jfi “d " Efl+i}- 


ii = R .- - R . . 
1 3.1 


2 R . 

^ 0 ^ R jC n 3R 


ij ^ ~ lj2, . .,,h~S 


f 3 hi - o .^ ♦ * » ^ 


,11 


-Ij = ( 1 >- 1 ) H ^ - S H + E , 


.n 


1% ®id ®h-13< i = h+3..., 

i (ll—l) R-i “* 2 R 4- T? 

^ 3 ^ %j» 3 ~ h+3,...jH 

^-Ij " ^ h + 13 ’ 3 R + 3 j...,n 


-3 ^h+2 3 » 3 - R+3,...,n 

"" i < 3, i,3 h+3,...,n 


^ h-lh - %- ih “ 


la-l 


n 


~ E+l ^ 1 - 


l-’’ 1 < 1 < d < n®iR 1 < , f ^ «ii 

^ 5 . ^ < 3 < n -^3 


h -2 


i=l i~il-l < 3 < k < / 3 l£“iEi^®ih-l+®ih+Sl)+(»- 2 )(Hi,_l+; 
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^st ^st "■ 


21 


2 

¥(¥ihrj 


fi.- 


R. . 


1 < i < 3 < n ^+1 < i < 3 < n >3 


h~2 


iil (Ei3 + Eit + El) - - H, 


h-1 


n 


+ l+I % 


i=l 


^h+lli+2 \+lh+2- 


jL— 1 -<—4 

It IS merely a roatlne worl!: to oheo:; that such transfonnat. 
ions make the matrix B to an identity matrix aiid f to 
Ihus if B-^ > 0 then l" is optimal. But since the choice of 
a(H) is not unique for the given point X", such condition is 

not unique, and hence the condition given is not necessary but 

sufficient. 


f s h— Ij hj t — h+1, ii+2 


-r 


EpSrr . ^ ^ £ R. 

i < 3- < 3 < n r+1 < i < j < n 13 


i-ills COIBpl-O'fce^ "fellB pjuOOf ♦ 

As discussed above, one can conceive of several other 
special cases for -.Yhich exact lower bounds can be obtained, 
vre, however, do not discuss all these oases here. 



Chapter Y 


LOWER BDUEL OP HIGHER DEGREE 


5 , 1 Introduction 

m 

So far we have restricted our scope of discussion to 
Bonferroni lower bound of degree two* In this concluding 
chapter, we give extensions of some results for the cases of 
lower bound of higher degree. In section 5.2, we generalise 
the class V* for the case of lower bound of degree v. Section 
5.3 deals \/ith a generalisation of the results obtained in 
chapter III, Since the concepts used in this chapter, are 
nothing but the generalization of the results given an the 
earlier chapters, here we do not discuss these results elaborately 


5. 2 Bonferroni laoer hound of degree v 

Ihe Bonferroni lower bound problem of degree v can be 
formulated as a linear program as follows; 

To maximize, 

II 

2 X. P. - E =^i i ^i i + 

i^=l H H 1 < i^ < ig < n ^12 

« ■ 

• • • + ( ~1 ) E j i ^ i i i 

l<i^<ig<,..<a^<n -^1^2 

(5.2.1) 
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sulDject to Z Xj - 2 

G J t ^ ig 

H.lg G J 


X.- + *.* 

tl^2 


... +(-l)''~^ Z i i 

1 < < ig < ... < < n 


< 1 


(5.2.2) 


for all J S IT. 


The abo-ve linear program IP (1) is of Z G.(= m^, say) 

i=l ^ 

variables with 2^-1 constraints, let P^C®) the set of 

feasible region of ItvCN)* -^s the previous case, we define 
a point X 6 B. with couponents x. , x. - ,...,x,. ^ - as 




..1, 


X (x^jXg, •• ‘ ***»^123- . .v, * ' * »^(n-v+l) . . .n 

We recall that the class Y* introduced in chapter II 




defines a subset of feasible points determined by integers h _> 1, 
a^, i Gh. Such a class of points can also be defined for the 
general case as follows, 

let us consider the polynomial 

n n n 

(h - Z aJCh- 2 a. + 1) ... (h ~ Z a- + V'-l) (5.2.3) 

i=l ^ i=l ^ i=l ^ 

where the variables are h, a-, i 6 1. This polynomial contains 


terms liXe t h^ , t hP a^^ a^^ ,»«.,t a-^ a.^. 

J-1 1 Xp X-1 Xc^ 


»3i 


Q.V 
iv » 


where, t, p-, q,:; are constants- 
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We define the class Y , JUce Y as, for a given integer 
vector (a,li) e h > 1 




^1 h( li+l ) ( n+2 ) « . « ( h+v *-l ) 


~ sum of the terms of the typ e 


t rf in (5.2.3):J, > o 

1 

X- ^ == — ' — — I sum of the terms of the type 

H^2 h(h+l)(h+2). ..(h+v-1) '■ 

_ Q.p 

t a^ in (5.2,3) > Q, q^ > o 

and similarly 


X . 


3-1 3-2 •••iv h(h+l)(h+2) .. .(h+v-1) 


£'sum of the terms of the type 


0 li Ip % 

t ii a- 3,- » « «aj in (5,2»3) 1, q^ ^ 0, 1 ^ i ^ 

"^2 ’ 1 ' — ^ 

In theorem 5.1, we shovv that the class is a subset of 
the feasible region of 11^(1). 

P- Pj.?z * 

Any X G Y c. R is feasible to the set of constraints 


(5.2.2). 


m. 


Pjiqof : Substituting the value of X 6 R ii terms of (a,h) 

in the constraint associated with any subset jCI, v/e get 


(h- Z a.) (h“ E 3-1+1) ••• (b— Z a.-hj-l) 
. - ~ G J ^ ^ ^ ® 

h(h-i-l)(h+2) ... (h-h)~l) 


'sHi; 


1 0 . 
(5.2.5) 



It is clear that for any integers h > 1, a^, i G the 
inequality (5.2.5) is true. Hence the constraint J is satisfied 
hy X e V*''. 

Shis coH^letes the proof. 


A special subset of can be defined as follows, 
For a giiren SC-IT, for integer h, f ^ h _< s-v, 


' _ v(v-l) _ _ _ 

Hh (h-K,-l)(h-H,-^) ^ 2 12 


and similarly. 


v(v-l) ...2 1 


h(h+l) ... (h+v-,1) 


; , ¥ i^ < i^ < 


< 


and x- = x- 


^ 1^2 


* ♦ ♦ — X 


X^lg.. .iy 


• • • iy ® ^ * 

= 0, otherwise. 


It is easy to see that for v = 2 such a liind of point 
give rise to class Vg defined in (1.5.2), Shis special class 
is useful to get a lower bound of degree v , C293, C473. 

Substituting the value of X in terms of (a,h) in (5.2,1), 

We see that the computation of Bonferroni lower bound of degree 
V can be viewed as a polynomial of degree v of integer variables* 
At present we do not discuss the computational aspect of solving 
such problems. 
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5. 3 Extr&ne paints of F-^(B) 

In this section we give some general results regarding 
the extreme points of 1^(1) » 

Any constraint of (5,2 ,2) can be written as 

d(J)® X < 1 

m 

where d(J) G E is defined as 


such that 

= 1 If e J 

= 0 if ^ J 


and 


d,. 


H^2 


= (-1)3-1 1 < 3 <v. 


m^ 

I'or such given vector d(J) 6 E v^e define a vector 
m. i 

d-(J) G R ^ as the first m.(= E ^0.) compcaaents of d(J)* like 
i 1 3 

the previous case, v\?e introduce the ranks of different orders. 

i'or a family of subsets of H, a(l)j the rank of order i, 

1 < i < V is def ined as the maximum number of linearly 

independent vectors in the set {dj^(J) J J G a(S)) , and 

this rank is denoted by r^(a(d)). 

In theorem 5,2, we give a characterisation of extreme 
points of first kind of 1^(1). 
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The or e m 5 >2 : , 

Tor a gi-ven subset U, let integers h > 1, a^i e 1\ S 
be chosen such that 

(i) 1 < h < s~v 

(ii) a a collection a(N\S) of NXs such that 

|a(l\S) ) = £ . . 

i=l 

r^(a(N\S)) = £ 

i=l 

(iii) 0 < 9(1) < s~v+l, ¥ I G a(N\S) where e(I) = h - £ a. • 

i e I ^ 

Then the point 1 of Y*^ defined by X = (a, h) such that a.^ = 
i e S' is an extreme point of Tv(W)« 

The proof of the theorem is omitted here because of the 
fact that the theorem can be proved following the same arguments 
as used in theorem 3*3* 

Remark : Here it is to be noted that like any point in 

Ext (1), any extreme poinit of first kind of Ty(N) satisfies v 
consecutive blocks of any particular column of R(S), which is 
independent of degree v , 

Using the concept of rank of order i, we proceed to 
generalise the idea of extreme points of different kind. 

let for a given Sc N, and given integers a^,h, 1 e l\S 
the subsets of I\S are partitioned as %, Xg as given in 
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section 3.4. Let us assume that r^(%) < E and 

i=l ^ 


exists ^(11X3) d n such that 


i=l ^ 

Then the point (a,h) g T*^ with a^, =1, i g s can 
extended to generate an extreme point which is called an 
points of second kind. 'Such concept is useful to define 
different kinds of extreme points for the LPy(H), 

We, hovTOirer, do not pursue it any further. 



there 


be 

extreme 
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Here we give q matrices for the 
in table 2.3, 


problems 


u is cus ss d. 


Pr jD^em 1 



0.49 5 

0.244 

0.244 

0.244 

0.118 

0.322“] 


0.244 

0 .495 

0.244 

0 .244 

0.118 

0.322 j 

s: 

0.244 

0.244 

0.495 

0 .244 

0.118 

0.322 ' 


0.244 

0.244 

i 

0.244 

0.495 

0.118 

0.322 j 


0.118 

0.118 

0.118 

0.118 

0.431 

0 .284 


„0 .322 

0 .322 

0.322 

0.322 

0 .284- 

0.591_;^ 



Pr oblem 2 




'"O .225 

0.Q_^ 

O^OBO 

0,090 

0.040 

0.020 


0*070 

0 .24o 

0 *09 5 

0,105 

0.035 

0.015 


0.080 

0.095 

0.245 

0*115 

0.030 

o.olo 


0.090 

0.105 

0.115 

0.250 

0 .025 

0.005 


0 .040 

0.035 

0.030 

0.025 

0.160 

o.olo 

. 

^ 0 *0 20 

0.015 

0.0 lo 

0.005 

0.0 lo 

0.0 75 









Problem 3 




"o .385 

0.210 

0.190 

0.210 

0 .215 

0f04o“ 


o 

• 

o 



^ * 




0.390 

0.180 

o.2io 

0.220 

0 ,045 


0 .190 

0.180 

0i360 

0.215 

0,215 

0.030 : 


0.210 

0.210 

0.215 

0.390 

0 .220 

0 .055 


0.215 

0 . 220 

0.215 

0.220 

0.385 

0*025 


,0.040 

0.045 

0.030 

0 #035 

0.025 

0.200 ' 
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Problem 4 


0.565 0.393 

0.393 0 « 56 3 

O.I60 0.368 
0.264 0.221 

0.340 0.231 

0.138 0.0 70 


0.162 0.264 
0.268 0.221 
0.384 O.lO? 
0.107 0.43 7 

0.179 0.255 

0.161 0.0 79 


0.340 

OD 

to 

• 

0 

0.231 

0.0 •79 

0.179 

0.161 

0 .255 

0.0 79 j 

0.522 

0.109 

0 .109 

0 .220 j 
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